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Abstract 

We design a theoretic tree-based functional representation of a class of Feynman-Kac 
particle distributions, including an extension of the Wick product formula to interacting 
particle systems. These weak expansions rely on an original combinatorial, and permu- 
tation group analysis of a special class of forests. They provide refined non asymptotic 
propagation of chaos type properties, as well as sharp L p -mean error bounds, and laws 
of large numbers for [/-statistics. Applications to particle interpretations of the top 
eigenvalues, and the ground states of Schrodinger semigroups are also discussed. 
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1 Introduction and formulation of the main results 
1.1 Introduction 

The field of Feynman-Kac semigroups, and their interacting particle models is one of the 
most active contact points between probability, theoretical chemistry, quantum physics, and 
engineering sciences, including rare event analysis, and advanced signal processing. Whereas 
it is clearly out of the scope of the present article to review these applications, some of them 
are briefly pointed out in the further development of this introduction. For a rather thorough 
discussion, the interested reader is recommended to consult the pair of books @J[S], and the 
references therein. 

The common aspect of all these applications is that they represent path integration of 
free Markov particle evolutions, weighted by some potential functions. The models we con- 
sider belong to the general class of branching and interacting particle systems. Particles 
explore the state space as a free Markov evolution; during their exploration particles with 
low potential are killed, while the ones with high potential value duplicate. From a more 
probabilistic point of view, these evolutionary genetic type particle models can be inter- 
preted as stochastic linearization, or as sophisticated interacting acceptation-rejection type 
simulation techniques. 

During the last two decades, the asymptotic analysis of these models has been developed 
in various directions, including propagation of chaos analysis, L p -mean error estimates, cen- 
tral limit type theorems, and large deviation principles. The purpose of the present work 
is to develop algebraic tree-based functional representations of particle block distributions, 
stripped of all analytical superstructure, and probabilistic ir relevancies. These weak expan- 
sions rely on an original combinatorial, and permutation group analysis on a special class of 
classical and colored forests that parametrize naturally the trajectories of interacting parti- 
cle systems. In order to modelize and compute the corresponding path integrals, we develop 
an original differential and combinatorial calculus, and take advantage of the natural per- 
mutation group action on particle configurations to handle the combinatorial complexity of 
the problem. Wreath products of permutation groups appear in the process, and new closed 
algebraic formulae for labelled forests enumeration are obtained. We also discuss Hilbert 
series techniques for counting the forests relevant to our analysis, that is with prescribed 
number of vertices at each level, or with given coalescence degrees. 

Let us mention that forests, colored graphs and their combinatorics have appeared re- 
cently in various fields such as in theoretical physics and gaussian matrix integral models (H| , 
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renormalization theory in high-energy physics or Runge-Kutta methods, two fields where the 
structure and complexity of perturbative expansions has required the development of new 
tools To the best of our knowledge, their introduction in the analysis of Feynman-Kac 

and interacting particle models is new. The Laurent type integral representations presented 
in this article do not only reveal the combinatorial essence of interacting particle distribu- 
tions. They also provide refined non asymptotic propagation of chaos type properties, as well 
as explicit formulae for L p -mean error bounds, and laws of large numbers for [/-statistics. 

The article is divided into four main parts, devoted respectively to the precise description 
of Feynman-Kac particle models and tree-based integral expansions, to the proof of these 
formulae, to the combinatorial and group-theoretic analysis, and to the extension to path 
space models with applications to propagation of chaos type properties . 

In section II .2| we describe a general class of Feynman-Kac semigroups associated with 
some free Markov evolution, and with a collection of potential functions. We motivate 
these abstract mathematical models with a brief discussion on a particle absorption model 
arising in the spectral analysis of Schrodinger operators. The mean field particle models are 
presented in section lL3l Section fl. 41 is concerned with the description of the corresponding 
particle block distributions. The main results of the article are housed in section IT31 

Section [2] is essentially concerned with the proof of the tree-based Laurent expansions 
presented in section IT31 

Section |3] is devoted to the combinatorial analysis of the class of forests involved in the 
tree-based Laurent expansions. The analysis includes some exotic planar forests, the branchs 
of which are allowed to cross, referred to in the following as jungles. The first subsection, 
section starts with a short review on the vocabulary of trees and forests. In section l3T2l 
and section 13.31 we provide an algebraic representation of forests and jungles in terms of 
sequences on maps. Sections 13. 4| 13.51 and 13.61 introduce the algebraic tools relevant to the 
combinatorial analysis. Permutation group actions on jungles are introduced. From the 
Feynman-Kac modelling point of view, they correspond to the action of the permutation 
group acting naturally on particle systems. This action is studied in detail. For example, 
a wreath product representation of the stabilizers is obtained. Closed inductive formulas 
for interacting particle models path enumeration are derived in section 13.51 Finally, in 
section |3~T1 we design an Hilbert series method for counting forests and the other famillies 
of relevant combinatorial and geometrical objects. As an application of these results on trees 
and forests, some explicit calculations of the first order terms of the expansions obtained in 
section |21 are given in section l3~%l and a generalized Wick formula for forests is obtained in 
section 13.91 

Section0]is essentially concerned with the extension of the previous analysis to Feynman- 
Kac particle models on path spaces. The corresponding tree-based expansions involve a new 
class of forests, with two distinguished colored vertices. The combinatorial analysis of these 
objects is slightly more involved, but essentially follows the same line of development as 
before. This study is summarized in section 14.11 and section 14.21 The next two sections, 
section 14. 3| and section 14. 4( are devoted respectively to unnormalized Feynman-Kac semi- 
groups on path spaces, and to the corresponding tree-based Laurent type expansions. In the 
final section, section l4~5*l we use the differential forest expansion machinery developed pre- 
viously in this article to derive precise propagations of chaos type properties of interacting 
particle models, including strong expansions of the particle block distributions with respect 
to Zolotarev type seminorms, sharp L p -mean error estimates, and a law of large numbers 
for [/-statistics associated to interacting particle models. 
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1.2 Feynman-Kac semigroups, and some application areas 

To avoid any state space restrictions, and any unnecessary topological assumptions, this 
article is concerned with abstract mathematical models in general measurable state spaces. 
The main advantage of this general set-up comes from the fact that it applies without further 
work to the analysis of path space Feynman-Kac models, and their genealogical tree based 
interpretations. 

Before we continue with the body of the introduction, we already present a series of 
definitions, the benefits of which are reaped in the subsequent sections. We let (E n ,£ n ) n >Q 
be a collection of measurable state spaces. We denote respectively by A4(E n ), V(E n ), and 
Bb(E n ), the set of all finite signed measures on (E n , £ n ), the convex subset of all probability 
measures, and the Banach space of all bounded and measurable functions / on E n , equipped 
with the uniform norm ||/|| = sup Xn£En \f(x n )\. We let /%(/„) = f En Vn(dx n ) f n (x n ), be 
the Lebesgue integral of a function f n G Bb(E n ), with respect to a measure fi n G A4(E n ). We 
equip M.(E n ) with the total variation norm ||/U n || tv = SU P feB b (E n )-.\\f\\<i lMn(/)l- We recall 
that a bounded integral operator Q n from E n ^\ into E n , such that, for any f n G Bb(E n ), 
the functions 

Qn(fn) ■ X n -\ G E n -l H-> Qn(fn)(x n -l) = Q 

are £ n _i-measurable, and bounded, generates a dual operator /i n -i >— ► fJ-n-iQn from M.(E n ^\) 
into M(E n ), and defined by (/%-iQ n )(/n) := /%-i(Qn(/n))- 

For a bounded integral operator Q\ from Eq into E\, and an operator Q2 from E\ 
into E2, we denote by Q1Q2 the composition operator from Eq into E2, defined for any 
[2 G Bb(E2) by (QiQ2)(/2) := Qi(<?2(/2))- The tensor power Q® 9 represents the bounded 
integral operator on En, defined for any F £ Bb(En) by 

Q®i(F)(xl...,x n )= [ [QnixldyD ... Qn(x n ,dy«)] F{y 1 w ...,y n ) 
JeI 

We consider a distribution r/o on Eq, a collection of Markov transitions M n (x n _i, <ix n ) from 
E n -\ into £/ n , and a collection of £" n -measurable, and bounded potential functions G n on the 
state spaces E n . To simplify the presentation, and avoid unnecessary technical discussion, 
we shall suppose that the potential functions are chosen such that 

< inf G n (x n ) < sup G n (x n ) < 00 (1.1) 

We associate to these objects the Feynman-Kac measures defined for any function f n G 
Bt(E n ) by the following formulae 

Vn(fn) = 7n(/n)/7n(l) with 7n (/ n ) = E[f n (X n ) U <k<n G k (X h )] (1.2) 

In H1.2I) . (X n ) n >Q represents a Markov chain, taking values in the state spaces (E n ) n >o, 
with initial distribution 770 ° n Eq, and elementary transitions M n from E n -\ into E n . The 
choice of non homogeneous state spaces E n is not innocent. In several application areas the 
underlying Markov model is a path-space Markov chain 

X n = (X' , . . . , X' n ) G E n = {Eq x . . . x ££) (1.3) 

The elementary prime variables X' n represent an elementary Markov chain taking values in 
some measurable spaces (E' n ,£ n ). For instance, in macro-molecular analysis, the elementary 
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variables X' n represents the monomers in a directed chain X n , with polymerization degree 
n. In this context, the potential functions G n reflect the intermolecular attraction or repul- 
sion forces between the monomers, and the corresponding Feynman-Kac model represents 
the distribution of the polymer configurations in a solvent. In signal processing, and more 
particularly in filtering problems, the chain X n represent the paths of a given signal, and 
the potential functions G n reflect the likelihood of their terminal values X' n , with respect 
to the sequence of observations delivered by the sensors. In this context, the corresponding 
Feynman-Kac model represents the conditional distribution of the path of the signal, given 
the sequence of observations. To motivate this article, let us also mention an important, 
and more classical physical interpretation of these Feynman-Kac flows. In theoretical par- 
ticle physics, these models represents the distribution of a particle evolving in an absorbing 
medium, with obstacles related to potential functions G n , taking values in [0, 1]. In this con- 
text, the particle X^ evolves according to two separate mechanisms. Firstly, it moves from 
a site x n -\ £ E n -i, to another x n £ E n according to elementary transitions M n (x n -\, dx n ). 
Then, it is absorbed with a probability 1 — G n (x n ), and placed in an auxiliary cemetery 
state X^ = c; otherwise it remains in the same site. If we let T be the random absorption 
time, then it is not difficult to check that 

7n(/n) = nfn(X c n ) 1 T >„] and r] n (f n ) = E[/ n p£) | T > n] (1.4) 

In time homogeneous settings (G n = G and M n = M), we have 

7n(l)=IP(T>n)~e- An 

The positive constant A is a measure of the strength and trapping effects of the obsta- 
cle. By a lemma of Varadhan's, A coincide with the logarithmic Lyapunov exponent Ao of 
the transition operator Q(x,dy) = G(x) M(x,dy). Whenever it exists, the corresponding 
eigenfunction h of Q represents the ground state of the operator Q. 

These Feynman-Kac models can also be regarded as a natural time discretization of 
continuous time particle absorption models associated with some Schrodinger operator. To 
be more precise, let us suppose that M(x,dy) = M A (x,dy) represent the transitions of a 
continuous time Markov process with infinitesimal generator L, during a short period of 
time A; in the sense that 

L(f)(x) = Jim I [M A (f)(x)-f(x]] 

for sufficiently regular functions / on E. Let us also suppose that the potential functions 
have the form G(x) = G A (x) = e~ v ^ A , for some non negative energy function V. In 
this case, it is readily checked that the infinitesimal generator of the particle absorption 
transitions 

Q A (x,d y ) = e -n*)A M A (X)dy) 

is given by the Schrodinger operator 

L(/)(x) - V(x)f(x) = hm 1 [Q A (f)(x) - f{x)] 

Under reasonably weak assumptions, the pair (A, h) coincide respectively with the top eigen- 
value, and the ground state of the Schrodinger operator L v = L — V . The estimation of 
these spectral quantities has recently received a lot of attention in modern numerical physics, 
and quantum chemistry. 
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A natural, and common key idea in most of the Monte-Carlo approximation models is 
to simulate a population of walkers or particles mimicking the evolution of a particle in 
an absorbing medium. These models can be interpreted as a natural mean field particle 
approximation of the distribution flow r/ n . We end this section with a precise description of 
these nonlinear models. By the Markov property and the multiplicative structure of (|1.2ft . 
it is easily checked that the flow (r/ n ) n >o satisfies the following equation 

7] n+1 = $ n +lO?n) (1.5) 

The transformations <& n +i : V(E n ) — ► V(E n+ i) are defined for any pair (rj n , f n+ i) £ 
(V(E n ) x B b (E n+1 )) as follows 

$ n+1 (j7 n )(/ n+ i) = !?IL-^±lj^±lL with Q n+1 (x n ,dx n+ i) = G n (x n ) x M n+X (x n , dx n+1 ) 
rj n {Q n+ i{l)) 

Notice that the evolution semigroup (Q Pt n)o< P <n of the unnormalized distribution flow is 
given by 

VO < p < n 7„ = 7pQp, n with Q Pt7l = Q Pt n-iQn (1-6) 

and the unnormalized Feynman-Kac measures 7„ satify "f n (fn) = ln-i{Qn{fn))- Equiva- 
lently, they can be expressed in terms of the flow (rj p )o< p < n with the formulae 

7n(/n) = Vn(fn) x 7n (l) with 7n (l) = JJ %{G p ) (1.7) 

0<p<n 

for any f n £ Bb(E n ). Returning to the physical interpretation, this product formula shows 
that the Lyapunov exponent A is related to the asymptotic mean logarithm value of the 
potential G n = G with respect to the measures r\ n 



- ^2 l°gMG)^-A 



n 

0< P <n 

Much more is true. In the case where the homogeneous transitions M n = M are reversible 
with respect to some measure /i, the ground state represents the limiting fixed point measure 
of the distribution flow r/ n ; in the sense that for any bounded measurable function /, we 
have that 

\ E V P (Gf)/r, p (G) ~ fM(hf)/(i(h) (1.9) 

0<p<n 

1.3 Mean field interacting particle models 

A natural mean field particle model associated with the nonlinear Feyman-Kac flow ()1.5j) 

is the -E^-valued Markov chain £n = (£,n N )i<i<N with elementary transitions defined for 
any G B b (E%) by 

N 

E{F n " { e ) )\e\)=*n(rn(e\)) ^ with -(ei)=^E^ ^ 

i=l 

In other terms, given the configuration ci-i a ^ ran k ( n ~ 1)) t ne particle system £A at 
rank n, consists of N independent and identically distributed random variables with common 
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distribution <& n Nj^ijY The initial configuration £q N ^ consists of N independent and 
identically distributed random variables with distribution rjQ. Although the dependency 
of the model on N is strong, due to the mean-field nature of the model, we abbreviate 
(£,n N )i<i<N to (C n )i<i<N when N is fixed, excepted when we want to emphasize explicitely 
the dependency on TV" of the model, e.g. as in the end of the present section. 
Notice that 

$n (m (eH)) (dx n ) = 1 f E G n-i(C-i) M n (e n -i,dx n )] 

Lj=l G n-H?n-l) \i=l J 

so that the particle model evolves as a genetic type model with proportional selections, and 
mutation transitions dictated by the pair of potential-transition (G n —i,M n ). In numerical 
physics, the selection transition is often called reconfiguration of the population of walkers. 
In the context of path space Markov chain model (|1.3|) . it is important to notice that the 
particle interpretation model is again a path particle model. During the selection stage, we 
select a path particle according to its fitness, and the mutation stage consists in extend- 
ing the selected path with an elementary transition dictated by the Markov transitions of 
the elementary prime chain. In this context, the particle model can be interpreted as a 
genealogical particle evolution model. 

The approximation measures (rj^ , 7^ ) associated with the pair of Feynman-Kac mea- 
sures (rj n , 7 n ) are defined by the empirical occupation measures 

^ = m (#0) 

and the unnormalized particle distributions 7^ are defined for any f n 6 Bb(E n ) by 

lZ(fn) = ^(fn) X 7^(1) With 7^(1) = II <^ ( L11 ) 

0<p<n 

For a rather complete asymptotic analysis of these measures we again refer the reader to [3] , 
and references therein. In particular, for any f n £ Bf,(E n ), we have the following almost 
convergence results 

Jim r/^ (/„) = r) n (f n ) and lim nfi(f n ) = 7n(/n) 

Therefore, the particle interpretations of the pair objects (A, h) introduced in (|1.8j) and (|1.9j) . 
are defined by simply with replacing r] n by the particle occupation measure 77^ . It is also 
well known that 7^ is an unbias approximation measure of 7 n , in the sense that 

E(T^(/n))=7«(/n) 
However, it turns out that for any n > 1 

E(^(f n ))=E(f n (e n N ))^Vn(fn) 

This means that the mean field particle interpretation model is not an exact sampling 
algorithm of the distributions rj n . In practice, it is clearly important to analyze these 
quantities, and more generally the distribution of particle blocks of any size q < N 

P; = Law(e,..,C W )eW) 
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1.4 Particle block distributions 

For any pair of integers (I, m) G (N*) 2 , we set [I] = {1, . . . , I}, and [l]^ the set of mappings 
a from [m] into [I]. By \a\, we denote the cardinality of the set a([m]), and for any 1 < p < I 
we set 

[Z]H = {aG[Z]H : \ a \ = p} 

The number m — \a\ will be refered to as the coalescence number of a. To simplify the 
presentation, when I > m, we also denote by (m, I) the set of all 

l\ 



(l-m)\ 



one to one mappings from [m] into [I], and by Qi = (I, I) the symmetric group of all permu- 
tations of [I]. 

From the pure mathematical point of view, the measures are better understood when 
they are connected with the g-tensor product occupation measures (resp. the restricted ex- 
tensor product occupation measures) on E n defined by 

a€[N]M y ' Q ae{q,N) 

Notice that these measures are, by construction, symmetry-invariant. That is, for any 
F G B b (E q n ) 

(F) = ((*%ym) and (^ f (F) = (^ f ((F) sym ) 

with the symmetrization operator 

Vg > 1 F G B b (E%) ^ (F) sym := 1 £ D a F 

In the above displayed definition, the operator D b stands for the Markov transition from E n 
into itself, associated with a mapping b G [q]^, and defined by 

D b (F)(x 1 n ,...,xl) = F(x% 1 >,...,#)) 

for any F G B b (En), and (x 1 , . . . , x 9 ) G F^. In particular, we may assume, without restric- 
tion, in our forthcoming computations on g-tensor product occupation measures, that F is 
a symmetric function. That is, for all s G G q , D S F = F. We write from now on, B s h vm {EV) 
for the set of all symmetric functions in B b (E n ). 

Notice also the following symmetry property, essential in view of all the forthcoming 
computations. Since, conditional to the are i.i.d., for any F G B b {En) and any 

a, b G (q, N} 2 , we have 

E[F(C {1 \---,Cn {q) )\ =E[F{& 1 \. . 
For instance, we have that 

E (W) 09 (F)) = E (F (£, . . . , a)) = < q (F) (1.12) 
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As in the case of elementary particle blocks of size one, the precise analysis of these q- 
tensor product measures is intimately related to the unnormalized pair particle measures 
((7^) 09 , (7n )® q ) defined for any F G B h (E q n ) by 

{l%) &q (F) = (^) Gq (F) x and (^(F) = (^)^(F) x ( 7j f (1))« 

We consider now the nonnegative measure Q^ q , indexed by the particle block sizes, on 
the product state spaces En, and defined for any F £ Bb(En) by 

®% q (F) := E(( 7 »r«(F)) 

In the further development of this article, we shall prove the following formula ()2.4|) 

This simple observation indicates that the distribution flow (Q^,j)n>o nas a linear evolution 
semigroup, for any particle block size q < N. Much more is true. As another consequence 
of the results presented in [I], the sequence of distributions Qn,g converges, as N tends to 
infinity, to the distribution of q non absorbed, and independent particles (Xn l )i<i< q evolving 
in the original absorbing medium. That is, we have that 

lim Q£,(F) = 1 ®«(F)=E(F(X%\...,X%1) l r i>„ ... 1 T9 >„) 

for any F £ Bb(En), and where T l stands for the random absorption time sequence of the 
chain (Xf' )fc>o, with 1 < i < q. 

This article is mainly concerned with explicit functional expansions of the deterministic 
measures F^ g £ V(E n ), and G M{En), with respect to the precision parameter N. 
These (Laurent type) expansions reflect the complete interaction structure of the particle 
model. Roughly speaking, the fc-th order terms of these integral representations represent 
the l/A fc -contributions of mean field particle scenarios with an interaction degree k (see 
Thru. ED - 

To describe more precisely these functional representations, it is convenient to introduce 
the abstract definition of the derivative of a sequence of finite signed measures. 

Definition 1.1. We let (Q n )n>i £ Ai(E)^ 1 , be a uniformly bounded sequence of signed 
measures on a measurable space (E,£), in the sense that sup N>1 \\Q N \\tv < oo. We suppose 
that Q N strongly converges to some measure O £ A4(E), as N f oo, in the sense that 

VF e B b (E) lim N {F) = 9{F) 

ATfoo 

The discrete derivative of the sequence (0 )n>i is the sequence of measures (d@ N )N>\ 
defined by 

d@ N := N [e N -&] 

We say that Q N is differentiable, if d® N is uniformly bounded, and if it strongly converges 
to some measure dQ £ Ai(E), as N f oo. The discrete derivative of a discrete derivative 
d@ N (of a differentiable sequence) is called the second discrete derivative and it is denoted 
by 

d 2 e N = n [de N - oe] 

The discrete derivative d® N of a differentiable sequence can itself be differentiable. In this 
situation, the derivative of the discrete derivative is called the second derivative and it is 
denoted by d 2 Q = d(d@), and so on. 
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At this point, it is convenient to make some comments. Firstly, the strong form of the 
weak convergence of measures considered above corresponds to the s-topology on the set of 
measures, defined as the coarsest topology for which all functionals fj, G M(E) \— > G M, 
(p G Bb(E), are continuous. For a more detailed discussion on the s-topology in the context 
of probability measures, we refer the reader to the article [2]. A sequence of real numbers can 
be interpreted as a sequence of measures with a constant value, so that the above definition 
applies to any sequence of real numbers. We also notice that the Euler's derivation operators 
d k satisfy the usual linearity properties. Finally, we notice that a sequence 0^ that is 
differentiable up to order {k + 1), has the following integral representation 



0</<fc 

with sup^ v>1 ||5 fc+1 Ar ||tv < °° > an d the convention <9°0 = 0, for / = 0. 
1.5 Forest based expansion formulas 

For transparency, we will often use terms associated with botanical, and genealogical trees. 
The forthcoming measures expansions will be therefore given in terms of signed measures 
indexed by collections of trees, colored forests and other analogous combinatorial and geo- 
metrical objects. More generally, most computations on particle measures will be rephrased 
in the course of the article in terms of enumerative problems on tree-like objects. In order 
to avoid lenghtening the present introduction, their definition and study are postponed to 
the next sections, and we present here our main results in a classical set-theoretic langage. 

Let A n ,q =def. ([g]^) n+1 be the set of (n + l)-sequences of mappings a = (a p )o<p<n from 
[q] into itself. 

[ q ]^[q}^---^[q] a ^[ q }^[ q ] 

Notice that we write in bold the symbols for sequences (of maps, integers...) such as a. 
We let A„ j(? be the nonnegative measure valued functional on A n ^ q defined by 



A n , g : a G A n , q ^ A* q = [jIq D ao Qf q D ai . . . Qt q D an ) G M(E q ) (1.13) 

The Markov operator Df, can be seen as a coalescent, or a selection type transition. In this 
interpretation, the population (i„ , . . . , i^') results from a selection of the individuals 
with labels in the set b( [q]). In addition, arguing as in Q1.4|) . the non negative integral tensor 
product operators Q® q can be seen as the overlapping of an absorption transition from 
(E n _i U {c}) q into itself, and an exploration transition from (E n _i U {c}) q into {E n U {c}) q . 
We then have that 

Af hq (F) = E {F(X^\ X^ q ) l T a,X> n • • • lTa,o n ) (1.14) 

In the above display, (A^' c>1 , . . . , Xn ,c,q ) represents the absorbed Markov chain on (E n U{c}) q 
with transitions Qn q D an , and initial distribution r/Q q D ao . 

Notice that the measures inherit a remarquable invariance property from their set- 

s in+rnrhire fine rm + nrnl left action of the grOUp C, 
, f 2 *™ 



theoretic definition. Namely, let us introduce the natural left action of the group +2 on 
An t g defined for all a G A n ^ q and all s = (sq, s n +i) G Gg +2 by 



s(a) := (s a s 1 \ sia 1 s 2 , SnQnS^) 
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Then, for any F £ B^ ym (En) we have 

A^(F) = A^)(F) (1.15) 

The identity wouldn't hold if F wasn't a symmetric function. However, as already mentioned 
in Sect. II. 4( this is not a serious restriction as far as the determination of the g-tensor 
occupation measures is concerned. The action of Gq +2 on A n> q induces a partition of A nA 
into orbits or, equivalently, a partition into equivalence classes, where a and b are equivalent, 
a ~ b, if and only if there exists s £ Q r q l+2 such that a = s(b). We write J~n,q f° r the 
set of orbits or equivalence classes. As we shall see later, T n ,q is nothing but (up to a 
canonical isomorphism) a remarquable subset of the classical set of forests. If a € A n ,q, the 
corresponding equivalence class is written a. We also write Stab(a) for the stabilizer of a 
in Qq +2 - According to the class formula, the number of elements in the equivalence class a, 
written #(a) or #(a) is given by: 

M -s _ wr +2 

m ' |5ta6(a)|' 

where we write \Stab(a)\ for the cardinal of Stab(a) and, more generally, \S\ for the cardinal 
of any set S. 

Explicit formulas for #(a) and for the various quantities associated to the action of G T q l+2 
on An q such as the number of orbits (and much more) will be given later, and form, from 
the combinatorial point of view, one of the cores of the article. As it will appear, due to 
its general form, our analysis paves the way for a systematical combinatorial treatment of 
all mean field path dependent approximations -e.g. without restrictions on the number of 
particles at each step of the approximation process. 

We are now nearly in a position to give our main results on Feynman-Kac particle 
models, but firstly, to simplify the presentation, need to introduce a pair of definitions. We 
start with some more or less traditional, and simplifying multi index notation. For any 
sequences of integers p = (pk)o<k<n, and 1 = (lk)o<k<ni we write p < 1 if and only if pf. < 
for all < k < n. We write |p| for (po + ... +p n ). Assuming now that p < 1, we use the 
multi-index notation 

n 

(!)p = n (^)pfc ' p! = n pki ' and ^(^p) = n s ^ fc,pfc ) 

0<k<n 0<k<n k=0 

where the s(lk,Pk) are Stirling numbers of the first kind. The difference (p — 1), and respec- 
tively the addition (p + 1) of two sequences is the sequence (pk — lk)o<k<n, and respectively 
(Pk + h)o<k<n- When no confusions can arise, we write N, and q, for the constant sequences 
(N)o<i< n , and (q)o<i< n - We also write 1, and respectively 0, for the sequence of unit inte- 
gers, and respectively null integers. The above definitions are extended to infinite sequence 
of integers p = (pk)k>o £ N N , with a finite numbers of strictly positive terms. Any function 
a : N i— > N on the set of integers into itself, is extended to integer sequences p S N N , by 
setting a(p) = (a(p k ))k>o- 

For a € A n , q , we write |a| the sequence ( | a-i | )o<i<n ■ The coalescence sequence of a is the 
sequence q — |a| of the coalescence numbers of the ajS. The coalescence degree of a is 

n 

coal{a) := (|q| - |a|) = ^(q - 

i=0 
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The subset of A n ^ q of sequences a such that |a| > q — p is written A n>q (p). Notice that, if 
a ~ b, |a| = |b|, so that the notions of coalescence degree and coalescence sequence go over 
to the set J- nq of forests. In particular, notation such as |a| or coal(a) is well-defined. In 
view of fla il. 151 for any choice (ao, a n ) of a representative of a forest f in T n ^ we set 

A^ g = A* g and |f| = |a| 

The subset of T n ,q associated to A njq (p) is written !Fn,q(p)- 

Theorem 1.2. The non negative measure represents the distribution of the coales- 
cent, and non absorbed Markov chain introduced in associated with a sequence 

A = (Ak)o<k<n of independent random mappings from [q] into itself, with common dis- 
tribution -ffq Ylae[q]W 7^y^ a ' That i- s > f or a M F £ B s ^ m {En), we have that 

®% q (F) = A* q (F) = E {F(X^\ X^ q ) l r A,i> n ■ • • l T ^>n) 
In addition, for any 1 < q < N , we have the polynomial expansion 



^n,q In > jyk 
l<k<(q-l)(n+l) 

with the collection of signed, and weak derivative measures <9 fc Q n ,<j given by the formula 

d k o = V V g(|f| ' q ~ r) #(f) A f 

° Un ' q ^ ^ (q), f | ^ 

We shall present several consequences of this theorem, including explicit descriptions of 
the first two order terms in the polynomial expansion, and a new extension of the Wick 
product formula to forests. 

The functional expansions of the distributions rely on the analysis of tensor product 
particle measures on path spaces. These models can be studied following the same line of 
arguments as the ones we used for the distributions Qn, q - As a parenthesis, we observe that 
these two measures coincide for constant potential functions 

(Vn > Vx„ e E n G n {x n ) = 1) (Vn > < N P^ q = Q% q ) 

To deal with general potential functions, we need to re- normalize the particle distributions 
7„ , and to extend the previous functional analysis to the path-space distribution of the 
particle measures defined for any integer sequence 1 = (4)o<fc<m by the following formula 

7 l,iV = {l N )&0 3 {7 N fh ^ ... ^ {l N r i n 

The analytic functional structure of these path space models relies on natural combinatorial 
techniques on colored rooted trees, and forests. Our result basically read as follows. 

Theorem 1.3. For any q, the sequence of probability measures (P^+i q )N>q is differ entiable 
up to any order. Besides, at any order N > q we have the closed formula 

r n+l,q-Vn+l+ ^ tfk or n+l,q ^ N \_(N-q)/2\ ° r n+l,q 

l<k<[{N-q)/2\ 

for some collection of signed derivative measures d k P n:q , whose values can be explicitly 
described in terms of a class of colored forests, with maximal coalescent degree k. 
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As before, we shall derive several consequences of these expansions, including an explicit 
description of the first order term, and a Wick product formula on colored forests. The latter 
provides sharp estimates for L^-mean errors between the particle occupation measures 77^ 
and the limiting Feynman-Kac measures rj n . Incidentally, combining this Wick product 
formula with the Borel-Cantelli lemma we obtain the almost sure convergence result 

lim (^ fi(F n , g ) = 

iV^oo 

for any bounded symmetric function F on En, with j E F(x n , . . . , x n )r] n (dx n ) = 0. This 
result is an extension of the law of large numbers for [/-statistics obtained by W. Hoeffding jjj 
for independent and identically distributed random variables to interacting particle models. 



2 Particle measures expansions on forests 
2.1 A preliminary stochastic tensor product formula 

The link between the two, usual and restricted, tensor products measures, {f]n)® q and 
(rjn ) &q relies in the end on a simple observation, that will appear to be fundamental for 
all our forthcoming computations. Throughout this section, integers N > q > 1 and a 
mesurable state space E are fixed once for all. 
Consider first the surjection 

vr : [q]® x (q,N) — ► [N]® 
(s, a) 1 — > as := ao s 

Lemma 2.1. Let b € [iV]^, then the cardinal \'K~ 1 ({b})\ o/7r _1 ({6}) only depends on the 
cardinal \b\ of the image ofb. It is given by 

\n-H{b})\ = (N-\b\) q _ lbl (q) lbl 

Proof. Observe indeed that, with the same notations as above, for b and a fixed, the equation 
b = as can have a solution only if the image of b is contained in the image of a. Let us 
consider such a a. Since a is injective, it has a left inverse a' for the composition of maps 
such that moreover a' a acts as the identity on [q]. The equation has therefore a unique 
solution s = a'b. Since the number n p of injections from [q] to [N] the image of which 
contains a fixed subset of [N] of cardinal p < q is 

the lemma follows. □ 

The lemma translates almost immediately into a functional relation between the two 
tensor product measures. Notice first the identity 

D a D b = D ab 

that holds for any pair of mappings (a, b) £ ([q]^) 2 - Let u,v be any linear combinations 

ctiOi, (resp. Yl Pibi) of elements of [qy®, or, in other terms, elements of the monoid 

■iei jeJ 

algebra of [q]^. We extend linearly the definition of D, so that: 

D u ■■= ^2 a i D ^ 

iEl 
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and 

D U D V = D uv = n, I,/),, /, 

ieijeJ 

We denote now by m(x) the empirical measure associated with an iV-uple x = (x l )i<i<N £ 

E N 



i=l 

For any integer q < N, we also consider the empirical measures on I? 9 defined by 

ae[JV][«] 

m ( x ) 0<? = tatt Z V ll) .-,* o(,) ) 

Although the notation is self-explanatory, notice that, in the sequel, we will write simply x Q 
for (x a ^\ . . . ,x a ( q ' ) ) and 5 x a for £(3.0(1) riI «(,)). 

Corollary 2.2. We Ziawe 

m(x)® g = m(x) 09 D L ^ 



T N = — V 



Proof. In view of the following identity, that holds for any function F G B b sym (E q ), with 
</ < V- 

S x aD b (F) = D b (F)(x a W,...,x a W) = F{x ab ^,...,x ah ^) = 5 xab (F) 
the corollary follows from lemma 0~T1 and the identity 

1 



(N) q 



(N) p (N - = 1 where N > q >p 



□ 



Corollary 2.3. The tensor product measure m(x)® 9 has a Laurent expansion: 



N k^d k L a 

0<k<q 



where 

0<fc<g 



and 



d k L q = Y, s(p,q-k) £ 

9"*<P<9 1 jP a£[g]M 
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Notice that, in the expansion of d k L q , the double sum involves only elements a £ [q]^ 
the coalescence number of which is less or equal to k. 

The corollary follows from corollary 12.21 and from the Stirling formula 

(N) p = Y, s( P ,k)N k 

l<k<p 

For q > 1, and recalling that s(q, q — 1) = — ( ^ )> we get 



^ - if E 



a L = — > a H — - > a 



S 2 L 



2 g! 

. 1 v-^ (q-l)(q-2) 1 1 
„ - (q,q-2 )- Y a ~~ -T \ E a+ T^- E 



Notice in particular that, since the composition map 

(g,g> X ( q ,N)^(q,N) 
is a surjection with all fibers of cardinal q\, we have 

m(x) &q D d o Lq = m(x) 09 

so that, asymptotically in N 

m(x)^~m(x) 09 

More precisely, we have, assuming from now on that x is generic (that is, Xi 7^ Xj if i 7^ j): 
Corollary 2.4. T/ie following formulas hold for the expansion of m(x) 8 " 3 

N ||m(x)®« - m(x)^|| to = 2 — > 9(9 - 1) 

? 

||m(x) 0<? L> afcL J| to = ^ |s(p,g- A;) I 

p=q-k 

Here, || stands fort the total variation, that is, for any linear operator Q on Bb(E q ), 

\\Q\\ tv := sup \Q(f)\ 

ll/l|oo=l 

Recall that the Stirling number of the second kind S(q,p) is the number of partitions of [q] 
into p non-empty subsets, so that 

\[q]$ ] \=S(q,p)(q) p 
N«= S(q,k)(N) k 

l<k<q 

The corollary follows then from our previous computations and from the following lemma. 
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Lemma 2.5. Let Q = Y u p ]P 5^, where the u p are arbitrary real coefficients, then 



we have 

\\Q\\tv= Yl Kl 

0<p<q 

Proof. Indeed, we have I |Q||to < Yl \ u p\ (N) p S(q,p) by direct inspection, since ||<5 x a|| to = 

0<p<q 

1 for all a G [q]^- Let us write sgn(k) for the sign of Uk, and let us introduce the function 
<f> £ B b {Et) defined by: 

4>{yi,-,y q )= Yl s 9 n ( k ) s \( yi ,..., yq ) 

0<k<q 

where X(yi,...,y q ) := \{y\, . . . ,y q }\ and S^fa y ^ stands for the Dirac function. Then, 

Q(.4>) = Yl \ u p\ (N) p S(q,p), and the lemma follows. □ 

0<p<g 

2.2 Laurent expansions 

The goal of this section is to derive a Laurent expansion of the measures Qn,q £ M(En), 
Qn, g ( F ) ■= H(ln)® 9 ( F )), F G &b{En), with respect to the population size parameter N. 
The following proposition is fundamental. Recall from Sect. 11.51 that we write J~ n ,q f° r the 



set of equivalence classes in An, q under the action of the permutation group Gq +2 



Proposition 2.6. For any integers q < N , any time parameter n £ N and any F G 

B^ m {En) we have the Laurent expansion 

E ^*^^q{F) 

with the mappings ^n,q introduced in \1. lty) . 

Proof. Indeed, combining the definition of the particle model with corollary 12.21 we first 
find that 



E (WT (F) | C\) = W(D) 9 x e (Vr (jo | e£3 

= (^(lll'xEl^f'^fF) 
= (^(l^xEktf.!)^^) 



AN) 
Sn-l 



Sn-1 



Using the fact that, conditionnaly to Cn^-i an( ^ ^ or ^ e B^ ym (En) 

a frA wrl = (eiT(ggg) = (7^T(ggg) 
ni " n " j 1 j (eo^(^(i)) (7^o^(Q^(i)) 

and 

(Wl)) = (7^-iQn(l)) 9 = (t£.i(<?„-i)) 9 = (T^l)) 9 
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we arrive at 



(t£-i)*W(i)) 



AN) 



AN) 

Sn-l 



= E ((7^-1) 89 (QT d l%f) 

Integrating over the past, this yields that 



AN) 
Sn-l 



(2.1) 



Using a simple induction, we readily obtain the formulae 



E 



E {{(^r q Qf q D L » . . . QTD l n) (F) 
E (Jjlo 9 F> L NQf 9 D L N . . . QTD L x) (F) 
1 V ^ A* (F) 

(n+l) ( q ), , ^n,q\ r ) 



aeA n 



a| 



where ^4 ni9 is the set of sequences of maps introduced in section 11.51 and where we use the 
multi- index notation introduced in the same section. In view of equation ()1.15j) . we know 
that 

A h n jF) = A^(F) 

whenever a and b belong to the same equivalence class in A. n ,q, and the Proposition follows. 



□ 



We are now in position to derive Theorem 11.21 



Theorem 2.7. For any time horizon n £ N, any population size N > 1, and any integer 
1 < q < N, we have the Laurent expansion 



&(*0 = & + E 



l<fc<(<?-l)(n+l) 



N k 



d k Qn, q { F ) 



(2.2) 



with the signed measures d k Q n ^ q such that 

d k Q n , q (F) = E E 

r<q, |r|=fc fe^n, 9 (r) 

for allF eB s b ym {E q n ). 

Proof. Indeed, according to Proposition 12.61 

OA - 1 V W 

J Arg(«+i) ( Q ) 

pe[l,g] n+1 VM; 



g(|f|,q-r) #(f) f 

(q)|fi 



Using the Stirling formula 



p fe^ n ,«,|f|=P 



(N) p = e •'KP' 1 )^' 11 

1<P 



E #(f) 



(2.3) 
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we also have that 

0*0 = E E s (p>!) ^itt^- E #( f ) A U F ) 

l<ql<P<q KH,P fe^„, g ,]f|=p 

from which we conclude that 

<0*) = E E ^-^^r E #( f )A^) 

r<qq-r<p<q Wp fe^ n>g |f|=p 

Observe then that #(f) = (q\) n+1 if |f| = q since, in that case, f is nothing but the 
equivalence class of all sequences of bijections in A n ,q. This ends the proof of the theorem. 

□ 

Corollary 2.8. Assuming that \\A^ \\t v < 1 for any f £ 3~n,q> the following formula holds 
for the expansion ofQ^ q 

lim N \\Q^ q -^% v <(n + l)q(q-l) 

Proof. We have (see proof of theorem I2.7|) that #(f) = (q\) n+1 if |f | = q and besides, 
|{r, |r| = 1}| = n + 1, and for any such r, s(q,q — r) = —(f). In conclusion, the case |f| = q 
contributes to (n + l)q(q — l)/2 to the asymptotic evaluation of N \ \Qn, q — ln q \\tv 

Now, let us consider a sequence r with |r| = 1, e.g. the one, written r^ with r, = 1, for a 
given < i < n. In that case, there is a unique f with |f| = q — r, which is the equivalence 
class of all sequences of maps a in A n ^ q with aj £ Q q , j / i and \ai\ = q — 1. In particular 
#(f) = (q\) n+1 ( q 2 ), so that, on the whole, the case |f| = Ti contributes to q{q — l)/2 to the 
evaluation. 

The corollary follows then from the Laurent expansion of Theorem 12.71 □ 

We conclude this section with noting that the above functional expansions also apply to 
the dot-tensor product measures 

More precisely, by definition of the particle model we have that 

^((^f q (F)) = E((e/ ? Qr(f))=QL, ? (« (2-4) 

3 Combinatorial methods for counting forests and jungles 

In the present section, we face the problem of computing the cardinals #(f) involved in 
our Laurent expansions. We also derive various identities relevant for the fine asymptotical 
analysis of mean field particle interacting models, such as the number of classes of elements 
in F n>q with a given coalescence degree. 

As already pointed out in the Introduction, we have choosen the langage of botanical and 
genealogical trees, both for technical reasons (since forests appear to be the most natural 
parametrization of elements in J~ n ,q)i an d also for the clarity of the exposition. 
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3.1 Some terminology on trees and forests 



In this section, we detail the vocabulary of trees, that will be of constant use later in the 
article. 

Definition 3.1. A tree (respectively a planar tree) is a (isomorphism class of) finite non- 
empty oriented connected (and respectively planar) graph T without loops such that any 
vertex of T has at most one outgoing edge. Paths are oriented from the vertices to the root. 
There is no restriction on the number of incoming edges. The empty graph is viewed as 
a tree and is called the empty tree. The set of all trees (resp. planar trees) is denoted by T 
(resp. VT). The vertices of a tree without incoming edges are called the leaves; the vertices 
with both incoming edges and an outgoing edge are called the internal vertices or the nodes; 
the (necessarily unique) vertex without outgoing edge is called the root. 

Definition 3.2. A forest f is a multiset of trees, that is a set of trees, with repetitions of 
the same tree allowed or, equivalently, an element of the commutative monoid (T) on T, 
with the empty graph Tq = as a unit. Since the algebraic notation is the most convenient, 
we write 

f = J 1 " 1 ! J^ m k 

for the forest with the tree Tj appearing with multiplicity mi, i < k. When T\ ^ Tj for i ^ j, 
we say that f is written in normal form. A planar forest f is an ordered sequence of planar 
trees. Planar forests can be represented by noncommutative monomials (or words) on the 
set of planar trees. The sets of forests and planar forests are written T andVT. 

For example, if T[ and T' 2 are two planar trees, f = T[T 2 T[ is the planar forest obtained 
by left-to-right concatenation of T[ , T 2 and another copy of T[ . 




Figure 1: a graphical representation of a planar forest f = T[T 2 T[. 

We say that a tree T is a subtree of T', and we write T C T', if the graph of T is a 
subgraph of the graph of T' , and if the root of T is also the root of T' . In a more pedantic 
way, a subtree of a tree is a connected subgraph of T' containing the root. A subforest f cf 
of a forest f is defined accordingly, as a collection of pairwise disjoint subtrees of the trees 
in f. 

The distance between two vertices in a tree is the minimal number of edges of a path 
joigning them. The height of a vertex is its distance to the root. We also say that a vertex 
with height k is a vertex at level k in the tree. The height ht(T) of a tree T is the maximal 
distance between a leaf and the root. The height ht(f) of a forest f is the maximal height 
ht(T) of the trees T C f in the forest 

ht(f) = sup ht(T) 

Tcf 
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For the unit empty tree, we also take the convention ht(Q) = — 1. 

We write ffe(f) for the number of vertices in a forest f at level k > 0. The number of 
vertices at each level in a forest is encoded in the mapping : 

v : f € v(f) = K(f))fc>o 
Notice that for any pair of forests (f , g) G J 7 , we have that 

v(ig) = v(f) +v(g) 

We write V for v(J-), which coincides with the set of connected initial integer sequences, 
that is integer sequences p = (pk)k>o £ ^ N satisfying the following property 

3ht(p) G N s.t. inf p^. > and sup Pk = 

k<ht(p) k>ht(p) 

For the unit empty tree To = 0, we use the conventions 

u(0) = = (0,0,...) and ht(0) = -1 

For any p G V, we write T p for the set of forests f such that v(f) = p. Our notation is 
consistent, since the height ht(f) of a forest f 6f p clearly coincides with the height ht(p) 
of the integer sequence p. For the null sequence 0, we have that Tq = {0}. 

When p G V is chosen so that po = 1, the set .Fp reduces to the set Tp of all trees T 
such that v(T) = p. The notation on trees, and forests goes over to planar trees, and planar 
forests in a self-explanatory way. For example, for any sequence p G V, VT P is the set of 
planar forests f, with ffc(f) = pk vertices at level k > 0. 

For any n G (NU {—1}), we denote by V n C V the subset of sequences p such that 
ht(p) = n. Notice that V_i = {0}. In this notation, the sets of all forests, and trees with 
height n are given by the sets 

Hn) ■= II Tin) := ]J T p 
P&Vn pev„,po = i 

The shift operator 

B : p= (Pk)k>0 G V ^ B(p) = (qi)i> G V, q t := p l+1 

induces a canonical bijection, still denoted by B, between the set of trees T p (resp. planar 
trees VT p ), with p G V„ + i and the set of forests Fb{p) (resp. planar forests VTb( p )), 
obtained by removing the root of the tree 

B : T G T p (resp. VT p ) ^ B(T) G F B{p) (resp. VF B[p) ) 

From the graph-theoretic point of view, the number of coalescences Cfc(f) at each level 
k > in a forest f G T is defined by the mapping 

c : f G c(f) = B(u(f)) - |f| 

The sequence c(f) is called the coalescence sequence of f. In the above displayed formula, 
the sequence |f| = (|f|fc)fc>o represents the number of vertices at level k > minus the 
number of leaves at the same level, that is the number of vertices with an ingoind edge, 
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whereas B(v(f)) k stands, by definition, for the number of edges between the levels k + 1 
and k. Finally observe that for any pair of forests (f , g) G J 7 , we have 

( |fg| = |f| + |g| and B{v({)+v(g))=B{v({)) + B(v(g)))^c(fg) = c{{) + c(g) 

We also say that the coalescence order of a vertex in a tree or a forest is its number of 
incoming edges minus 1. The coalescence degree of a tree or a forest f is the sum of the 
Cfe(f) or, equivalently, the sum of the coalescence orders of its vertices. We say that a tree 
is trivial if its coalescence degree is 0. 

3.2 On the algebraic structure of trees and forests 

In a planar forest f , vertices at the same level k > 0, are naturally ordered from left to 
right, and therefore in bijection with [^(f)]. Planar forests (resp. planar trees) f G VT 
of height ht(f') = (n + 1) are therefore canonically in bijection with sequences of maps 
a = (ao, ...,a n ), where a k is a weakly increasing map from \pk+i] to \p k ], with k = 0, ...,n, 
and p G V n+ i (resp. and po = 1). 

To clarify the presentation, we write a f to indicate that a given mapping a between 
sets of integers is weakly increasing. In this notation, for any p G V n+ i we have that 

VT P ~ IA P : = {(a , ...,a n ) G ([p ] [pi] x . . . x \p n ]^+^ : V0 < k < n a k |} 

We write PP(a) (resp. PT(a)) for the planar forest associated to a sequence in IA P (resp. 
with po = 1). In this interpretation, when f G VT P is represented by a sequence of 
mappings a = (ao, a n ) G X»4 p , the number of coalescent vertices at level k is given by 

Cfc(f') = v k+1 (f') - \a k \ 

Notice that the definition is coherent with the definition of the coalescence sequence associ- 
ated to a. In addition, when the integer sequence p G V n+ i is chosen such that po = 1, the 
mapping B introduced in the previous Section can be rewritten as follows: 

B : PT(a) G VT p ^ 5(PT(a)) := PF(5(a)) G VT B{p) 

with the shift operator P(ao, • • • , a n ) := {a\, . . . , a n ). These definitions and notations can 
be extended to arbitrary sequences of maps and arbitrary trees and forests as follows. 

Definition 3.3. Let a = (ao, ai, a n ) be any sequence of maps a k : [pk+i] — ► \pk\> w ith 
k = 0, . . . ,n. To a sequence a is naturally associated a forest P(a): the one with one vertex 

for each element o/LKiofafc]* anc ^ a ec ^ e -^ or eac ^ ^ r (*' a fc(*))'^ e [Pfc+i]- ^ e sequence 
can also be represented graphically uniquely by a planar graph J(a), where however the edges 
between vertices at level k + 1 and k are allowed to cross. We call such a planar graph, where 
paths between vertices are entangled, a jungle. The set of jungles is written J . 
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Figure 2: The entangled graph representation of a jungle with the same underlying graph 
as the planar forest in Fig. 1. 

Notice that planar forests are particular examples of jungles, so that, for any sequence 
of weakly increasing maps a we have PF(a) = J (a). Notation such as T p , with p G V n+ i, 
extends from forests to jungles in a self-explanatory way. Jungles j G J p of height (n + 1) 
are therefore canonically in bijection with sequences of maps (ao, ...,a n ), where a& is a map 
from [pfc+i] to [pk], with k = 0, ...,n. More formally, we have that 



J P ^A P := ([po] N x...xW M ) 



To make these algebraic representations of planar forests and jungles more precise, it is 
convenient to introduce another round of notation. For any integers p,q,r,s > 1, and any 
pair of maps (a, (3) G ([q]^ x [s]^), we let (a V/3) be the weakly increasing map from [p + r] 
into [q + s] defined by 



(«V/3)(») 



a(i) if i € [p] 
q + (3(i) if % = p + j with j G [r] 



If a and /3 are weakly increasing, so is a V (3. For null integers p = 0, or p = = q we use the 
conventions [g]^ = {0}, and the composition rules (0V0) =0, (a V0) = a = (0Va). Observe 
that for any (cr,r) G (£ ? x Q s ), any (a,/3) G {[q] [p] x [s]M) and any (p,9) G (£ p x C? r ) we 
have the composition formula 

(a V r) o (a V (3) o (p V 0) = (a o a o p) V (r o (3 o 0) (3.1) 

For any p G £? p + r with a and /3 as above, we also quote the following factorization rule 

(a V f3) o p G V =>■ p G (<? p V Gr) 

To check this final assertion, we suppose there exists some i G [p] such that p{%) = p + j, for 
some j G [r]. In this case, we would find the following contradiction 

(a V (3){p{i)) = (a V (3)(p + j) = q + f3{j) [q] 

By a simple induction argument, for any «j G = 1, ...,m, we also have that 

(cui V ... V a m ) o p G ([gi] [pi] V ... V [g m ] W ) => p G (£ P1 V . . . V S P J (3.2) 

for any permutation p G ^ Pl +...+ Pm , and m > 1. 

More generally, for any sequences of integers p G V, and r G V, and any sequences of 
maps a = (ak)k>o £ A p , and (3 = ((3k)k>o £ A r we denote by (a V (3) the sequence of maps 

a v /3 := (a k V /3 fc )fc>o G -4 p+r 

Notice that, if a G TA P and /3 G X4 q , then a V /3 G Z„4 p+q . In the above displayed formula, 
we have used the conventions = 0, and (3i = 0, for any k > ht(p), and I > ht{r). 

Recall that sequences of weakly increasing maps a G IA P , and (3 G IA r can be rep- 
resented graphically by the planar forests PF(a) G VT P , and PF([3) G VT T . In this 
interpretation, the composition formula on the noncommutative monoid VT is now given 
by 

PF{a)PF{[3) = PF(a V 0) G 7^ p+r (3.3) 
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The same constructions hold for non necessarily increasing maps. In this situation, the 
sequences of maps a G A p , and (3 G A r can be represented graphically by the jungles 
J (a) G J p , and J(f3) G J r: and we have the (noncommutative) product formula on the set 
J 

J(a)J(/3) = J(a V /?) € J p+r (3.4) 
3.3 Automorphism groups on jungles 

For any given sequence of integers p 6 V n+ i, the product permutation group 

G P ■= (G P0 x G Pl x ... x G Pn+1 ) 

acts naturally on sequences of maps a = (ao, a±, a n ) G A p , and on jungles J(a) in J p by 
permutation of the vertices at each level. More formally, for any s = (sq, G G P this 

pair of actions is given by 

s(a) := (soao^r 1 , •••, s na„s~| 1 ) and sJ(a) := J(s(a)) 

An automorphism s G G P of a given jungle J(a) G J p is a sequence of permutations that 
preserves the jungle, in the sense that sJ(a) = J(a). The set of automorphisms of a given 
jungle J(a) G J p coincides with the stabilizer of J(a) with respect to the group action. 

Stab(J(&)) = {s G Q p : sJ(a) = J(a)} 

By the definition of a jungle, we also have that 

J(a) = J(b) a = b 

from which we find Stab(J(a)) coincide with the stabilizer of a with respect to the group 
action. The stabilizer theorem also provides a one to one correspondence ip a between the 
orbit set 

Orb(J(a)) := {sJ(a) : s G G P } C J p 

and the quotient G P /Stab( J(a)) 

V? a : J(b) G Or6(J(a)) ^ s ( ^ h) Stab(J(a)) G g p /Sta&(J(a)) 

for any choice of S( a>b ) G ^ p such that S( a b ^(a) = b. 

Notice that if two sequences a and b differ only by the order of the vertices in J(a) and 
J(b), that is by the action of an element of G p , then the associated forests are identical: 
F(a) = F(h). Moreover, the converse is true: if F(a) = F(b), then J(a) and J(b) differ 
only by the ordering of the vertices, since they have the same underlying non planar graph. 
In this situation, a and b belong to the same orbit under the action of G p . In particular, 
we get a result already mentioned in the Introduction 

Lemma 3.4. The set of equivalence classes of jungles in Jp under the action of the per- 
mutation groups Qp is in bijection with the set of forests Tp. In particular, due to the 
equivalence between the two notions of jungles and sequences of maps, it follows that the set 
^n,q of equivalence classes in A n ^ q under the natural action of Gq is canonically in bijection 
with the set of forests with q vertices at each level < k < n + 1. 

Here, q denotes, once again, the g-constant sequence of lenght n + 1. 
Since the root mapping ao of a planar tree is constant and equal to 1, we finally observe 
that for any a = (ao, . . . , a n ) G IA P , we have 

Stab{J{&) = Stab(PT(&)) = {U} x Stab(B(PT(a))) 
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3.4 Orbit sets of planar forests 

From the Feynman-Kac mean-field approximation point of view, jungles in A n ,q describe 
trajectories of families of particles with prescribed coalescence properties. In particular, two 
trajectories that are equivalent under the action of have the same statistical properties. 
In the present section, we face the general problem of computing the cardinals #(f), defined 
as the number of jungles in f G F p , for some integer sequence p G V n+ i, when f is viewed 
as an equivalence class according to the previous lemma. Let us write f in normal form 
(with the TjS all distincts) as a commutative monomial of trees: 

f = T 1 mi ...T™ fc , mi + ... + mk =p 

Let us also choose arbitrary representations T( 6 VT of the TjS as planar trees, so that we 
can view f as the forest associated to the planar forest 

f = (T[) m \..(T' k ) mk 

obtained by left-to-right concatenation of m\ copies of T[, 777-2 copies of T^,..., copies of 
T' k . We write T" --T^ for the expansion of the noncommutative monomial 

f = {T[) mi ...{T' k ) m " =T{'...T^ Q 

as a product of planar trees (without exponents). 
From previous considerations, we have that 

#(f) = #Orb(f) = # {g € J p : 3s £ G P s.t. g = sf } 

Due to the class formula, we also know that 



#(f) = #Orb(?) 



\G P \ p! 



Stab(i>) Stab(f') 



where we remind that Stab(f') stands for the stabilizer of the jungle f 6 J p . 

The computation of #(f) will be done by induction on (n + 1) := ht(f). Let us assume 
that we know #(g) for any forest of height less or equal n. Notice that, then, we also know 
#(t) for any tree t of height (n + 1), due to the canonical bijection B between trees and 
forests: #(t) = #(Bt). From the previous discussion, the problem amounts to compute the 
cardinals of the stabilizers Stab(f') inductively with respect to the height parameter. The 
following technical lemma is instrumental. 

Lemma 3.5. There is a natural isomorphism 

Stab(T{'...T^) n ({1 P0 } x G {pu ..., Pn+l) ) ~ Stab(T{') V ... V Stab(T^) 

where l po stands for the identity in Q po . 

Proof. We let a 1 = (a l , . . . , a l n ) be a sequence of weakly increasing maps such that T" = 
PT(a l ), with I <i <pq. In this notation, we have that 

f = T['...T^ = PTia 1 ) . . . PT(a p ») = PF(a) 

where PF(a) stands for the planar forest associated with the sequence of weakly increasing 
maps 



a = (a 1 V . . . Va w ) = (a ,. . . ,a n ) with VO < I < n a t = a] V . . . V a\ 



pa 
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Notice that the stabilizer of the planar forest PF(a) is given by 

Stab{PF(a)) = {s G G p : s(a) = a} 

and for any s = (sq, . . . , s n+ i) G Q p we have that 

s(a) = a ^ (VO < I < n s t o {a} V . . . V af°) o = a/ V ... V of ) 

Therefore, using the composition formula (j3.1|) and the factorization rule following the 
formula, it is not difficult to check that any sequence 

s = (l po ,si, . . . ,s n+ i) G Stab(f) 

can be represented in an unique way s = (s^ V ... V s( p °)), in terms of some permutation 
sequences 

s» = (l 1)S » ) ... )S « 1 ) 

between the vertices at each level of the planar tree PT(a l ), with 1 < % < pq. In terms of 
the level parameter 1 < I < (n + 1), we also notice that 

Sl = V ... V S [ Po) 
To take the final step, we observe that 

s iai s T+ \ = (,« V ... V s^) o (aj V . . . V af°) o ((s^ V ... V (s^?)- 1 ) 
= ( S W o aj o V ... V (,<*»> o a f° o 

This implies that 

s(a) = (s« V ... V s^) (a 1 V ... V a Po ) = (s^a 1 ) V ... V s (po) (a Po )) 

from which we find that sJ(a) = J(a) if, and only if, we have 

^^(a 1 ))... J(s( p °)(a P0 )) = ^a 1 )... J(a P0 ) 
This is clearly equivalent to the fact that 

J(s«(a i )) = S Wj(a i ) = J(a*) 
for any 1 < i < Po- This yields that for any 1 < j < p^, we have 

s« = (l 1 ,4 i) ,..., S i i } 1 )G5ta6(Pr(a t )) 
from which we conclude that the mapping 

s G Stab(PT^) . . . PT(a™)) n ({l po } x <? {pi ,..., Pn+l) ) 

i 

( S W V ... V s^°)) G Stab{T(a})) V ... V Sia&(T( a P°)) 
is an isomorphism. This ends the proof of the lemma. ■ 
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The first coordinate mapping 

7T : s=(* ,...,Sn+i) e Stab({T{) mi ...(T^ m ") ^ir(s) = s G (<? mi V . . . V G nik ) 

is a surjective map (it even has a natural section, the construction of which is omitted). We 
therefore have that 

Ker{n) = Tr-^lp,,) ~ Stab(T{') V ... V Stab^J ~ Stab(T[) mi x . . . x Stab{T' k ) mk 
This yields the isomorphism formula 

5ta6((^) mi ...(^) m *)/ (Stab(T[) mi x . . . x Stab(T' k ) m *) ~ (0 mi V ... V £ m J 
from which we readily deduce the following recursive formula for the computation of #(f). 
Proposition 3.6. VFe /tawe 

fe 

l^a^^r 1 ...^)^)! = JJ (m*! \Stab {B{T>))D 

i=i 

3.5 An inductive method for counting jungles 

We conclude this section by giving a closed formula for #(f). Let us introduce for that 
purpose some further notations. 

Definition 3.7. Let f be a forest written in normal form f = T™ 1 ...T™ fc . In this situation, 
we say that the unordered k-uplet (mi,...,mfc) is the symmetry multiset of the tree T = 
B~ l (i) and write 

s(T) = s(B-\T™K..T^)) = (mi,...,m fc ) 

The symmetry multiset of the forest s(f) = T™ 1 ...T k mk is the disjoint union of the symmetry 
multisets of its trees 

( \ 



sCTf 1 ...^) 



s(Ti),...,s(Ti),...,s(r fc ),...,s(r fc ) 



\ mi terms m k terms / 

For example, the symmetry multiset of the first tree of the forest f displayed in Fig. 
1 is (1,1); the symmetry multiset of B _1 (f) is (2,1); and the symmetry multiset of f is 
(1,1,3,1,1). 

We also extend the map B from trees to forests to a map between forests and, if f = 

jmi _ T m k isaforest; get 



5(f) = B(T 1 ) mi ...B(r fc ; 



m k 



From the point of view of graphs, the operations amounts to removing all the roots and all 
the edges that have the root as terminal vertex from the graph defining f . The definition 
extends naturally to planar trees and planar forests. In particular, if f is a planar forest 
with f as its underlying forest, we set s(f') := s(f). 

Theorem 3.8. The number of jungles (or, equivalently, of Feynman-Kac mean-field type 
trajectories with the same statistics) in f G T v , with p G V n +i, is given by 

#(f) = -s— ^ 

n s(B*(f))! 

?=-i 

where we use the usual multi-index notation to define s(B l (f))\. 
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Proof. It clearly suffices to prove that 

n 

\Stab(f')\ = J] s(^(f'))! 



i=-l 



for any planar forest f of height (n + 1). We check this assertion by induction on the height 
parameter. First, we observe that a planar forest f of height 1 can be represented as a 
noncommutative monomial 

f = {T[) m \..(T' k ) m * 

with different planar trees T[ of height 1, and some sequence of integers m^. In that case, 
the planar forests B{T[) reduce to elementary planar forests with null height, with S(T!) 
elementary planar trees with null height. This yields that 

ht(T l ) = l^\Stab{B{Tl))\=s{T>)\ 

By proposition 13.61 we conclude that 

k 

\Stab(f')\ = J] ( mj \ (s(T>)\D = s(B-\tyl 8(f)! 

i=i 

This ends the proof of the formula at rank 1. Suppose now that the assertion is satis- 
fied for any planar forest with height n. By proposition 13.61 for any planar forest f = 
(T{) mi ...(T/) m * , with height (n + 1), and written in terms of distinct planar trees Tj, we 
have 



\Stab(?)\ = S (B-\f))\ H \Stab(B(T!))\> 



i=i 



Since the planar forests B(T() have height n, the induction hypothesis implies that 

n— 1 n 

\Stab(B(Tl))\ = H s(Bi(B(T>))) = Yls(Bi(T>))! 

j=-l j=0 



This yields that 



k / n \ n / k \ m 



\Stab(f')\ =s(B- 1 (f '))!]] n s (^'(^'))! =s(B- 1 (f'))!lI II s ( Bi ^) 

i=l \j=0 J j=Q \i=l / 

Recalling that B*(f') = B^(T{) mi ...B^{T' k ) m \ we also find that 

s(B^f')) = (s(B^T{)), S (B\T{)), . . . , *(Bi{li)), . . . , b(^(3^))) 
mi terms m fc terms 

We conclude that 

|5ta6(f / )|=s( J B- 1 (f'))!ns(^(f'))! 

i=o 

This shows that the formula is also true for planar forests with height (n + 1). The inductive 
proof of the theorem is now completed. ■ 
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3.6 Wreath product representation 

This section is concerned with an original wreath product interpretation of the stabilizer 
of planar forests. To be self-contained, let us recall (see e.g. [8]) that the wreath product, 
written Qk < G > of an arbitrary group G with the symmetric group of order k, Qk, is the 
group that identifies, as a set, with Qk X G \ and where the product is defined by: 

(s,gi, ...,gk) 0',3i,-,5fc) := (ss', ffi5s-i(i), -, gug's-nk)) 

In addition, if a group G acts on a set S, and if Qk acts by permutation on S , that is, if 
for any s E Qk and any (x\, —,Xk) E S* the action of s on (xi, ...,£&) is given by 

= (a; s -i (!),..., a? s -i (fc )) 

we then have, for any (<?i, g7c) G 

s(gi, ...,5fc)s _1 (xi, ...,x fc ) = s(pi, ...,gk)(x s (i), -,x s (k)) 

= s(gix s (i),-",gkx s ( k )) 
= (p s -i(i)^ij •••,5 f s -i(fc)2;fc) 

We are now in position to state and prove the main result of this section. 

Theorem 3.9. For any planar forest (T{) mi ...(TL) mk E VT V , with height (n+1) written in 
terms of distinct planar trees T[, for some sequence of integers p E V n +i, there is a natural 
isomorphism 

Stab{(T[) m \..{T' k ) m *) ~ [Q mi < Stab(T[) >) V . . . V (£ mfc < Stab{T^) >) 

so i/iai, in particular, the stabilizer of any jungle in the equivalence class of f is conjugated 
m Q p to (Q mi < Stab{T[) >) V . . . V (Q mk < Stab^) >). 

Proof. We write T"...Tp Q for the expansion of the non commutative monomial (T[) mi ...{T' k ) mk 
as a product of planar trees (without exponents). Let us also choose an arbitrary 

s = (s , ...,s n+ i) E Stab(?) 

This implies in particular that the root permutation sq E £/ po is of the form 

s = (4 1} v ... v 4 fe) ) e (S mi v ... v £ m J 

Also observe that so acts on the planar forest T"...T'' by permutation of the planar trees: 

S „[Tf...^J := (T': ;Htr .rU M ) 

This action by permutation can be viewed as a permutation (by blocks, according to the 
permutation of the planar trees) of the vertices of the planar forest, and this action can 
therefore be rewritten uniquely under this process as the action of an element 

so := (so, s ,i, • • • , S ,n+l) E Qp 
on the vertices of the planar graph of T"...T" 

s [T{'...T"] = so(T 1 "...T' / ) 
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Notice that s depends on the planar forest T"...Tp Q , since it acts non trivially on the vertices 
at levels i > 1. In general, the embedding so G Q po "-^ so G G p is different from the canonical 
embedding s G G po ^ (s , l Pn+1 ) G £ p . 

We consider then the map (of sets) 

tx : ste&(f) ^ (g mi v . . . v Grn k ) x (5tob(f 7 ) n ({i po } x s (Wr .., Pn+l) )) 

s = (s ,...,s n+ i) i ^ /i(s) = (s ,ss x ) 

In lemma 1331 we have proved that (Stab(T" ...T") C\ ({l po } x G( Pl ,... )Pn+1 ))) is isomorphic 
to the permutation group Stab(T") V ... V Stab(Tp Q ). Therefore, the permutation sequence 

SS = S 1 S 0,1> • • • » S n+l s 0,n+l) 

can be alternatively represented as 

ss 1 = (V 1 ) V ... V s( p °)) G Stab{T[') V ... V Stab(T^) 
Using these representations, we have that 

p( s ) = (( s « v . . . v 4 fe) ), (s (1) V...VsW)) 

= (4 1} , (s« V ... V s^))) V (4 2) , (s( mi+1 ) V ... V s ( mi+m2 ))) 

\y V (^Sq^ / V" 11 + '" +mfc ~ 1+1 ) V y s K+...+m t _ 1 +m t )jJ 

On the other hand, we have 

//(St) = (s to,St(s to) _1 ) 

= (s t , sttg Sq : ) = (s io, (ssq )so(ttg x )so x ) 
and for any level index 1 < I < u we find that 

scyfe^Ki = «o,/(*i (1) V...Vt{ po) ) So j 

= (^ 1)rl(i)) v ...vt; (41>rl(mi)) ) 

v ^((4 2) )^ 1 (™i+i)) v vt ((4 3 V 1 ("»i+"w))) 



This implies that 



v v j t ((4 ') ("'i+-+'"t-i+l)) v v ^(s^) ( m i + ---+ m fc-i+ m fc))^ 

s (ttg ^Sg 1 = (l po ,ri, . . . ,r n+ i) 



with a sequence of mappings 

(l P0 ,n, . . . , r n+1 ) = (r« V ... V r«) G Sta&(I?) V ... V 5to6(2^) 

given by 

r (i)v...vrW = t« s o 1) r 1 ( 1 ))v...vt« s o 1) )" 1 ( m O) 

r (mi+l) v ... V r ( mi+m2 ) = t(( s< o ) )' 1 ( m ^+ 1 )) V ... V t^ s o 2) ) _1 ( mi + m2 )) 
r (mi+...+m fc _ 1 +l) \/ ... V r ( m i+---+ m fe) = t^ s o fc) )~ 1 ( m i+---+ m fc-i+ 1 )) V ... V t^ s o fc) ) _1 ( mi+ --- +mfc ^ 
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This yields that 

(ss X ) o (s (tt x )so r ) = (s« V... Vs( p «') o (r« V... VrW 

= ((s (1) o r (1) ) V ... V (s (po) o r (po) )) 
We finally obtain the following formula 

Kst) = (4 1} o ((s« o t ((4 1) )- 1 W)) v ... v ( B ("*> o t ((4 1) )- 1 (™ 1 )) ) )) 

v (4 2) o t®, (V mi+1) o t ((4 2) )- 1 (™i+i)) ) v ... v ( s (™i+^ o t ((4 2) )- 1 (-i+™ 2 )))^ 

V... 

= { (4 1} , (s« V ... V S (™))) (4 1} , (t« V ... V t^))) } 

v { (4 2) , (s< mi+1 ) v ... v s ( mi+m2 )) ) (4 2) , (V mi+1 ) v ... v t ( mi+m2 )) ) } 



V... 

and the Theorem follows. 



3.7 Hilbert series method for forests enumeration 

In the present subsection, we face the problem of computing other cardinals relevant to the 
analysis of Feynman-Kac mean field particle models, according to our Theorem 11.21 For 
example, we want to be able to compute the number of forests in any ^" p ; the cardinals 
|^*n,g(r)h or the number of forests in T nA with a given coalescence degree. Notice, by the 
way, that the notions of coalescence sequence and degree, as defined in the Introduction for 
sequences of maps in A n ,q, and forests in T n ^ go over to arbitrary forests and trees. 
We let (5) be the commutative monoid generated by a subset of forests 5 € J 

(S) = {f™ 1 . . . f™ k : k > 0, and VI < i < k mi > and /< G S} 

with the convention /J 711 . . . f™ k = 0, for k = 0. Notice that there is a canonical map from 
(S) to J 7 , which is an embedding, e.g. if S is a set of (distinct) trees. We also consider 
a set of variables x(T) indexed by the set of trees T G T, and we let M[T] be the set of 
monomials defined as follows 

M[T] = {x(Ti) mi . . . x(T k ) mk ■ k > 0, and VI < i < k m, > and Tj G T} 

For any mapping x ° n the set of forests 

X : f G T X(f) G M[T] such that V(f , g) G T 2 X (fg) = *(f) *(g) 

we have the Hilbert series expansions 

- ^ x(g) n i _\(f\ = J2 x ^ ( 3 - 5 ) 



1- v(f) ^ AW 11 1- v(f) 

AV ; ge{f> fes AV y ge{5> 

We let V* be the subset of all multi-indices p G V, with height ht(p) < n 



VI : = U_i< fe<n V 



A.- 



30 



Definition 3.10. The (multidegree) Hilbert series of forests, HJ-(x), is the Hilbert series 
associated to the partition of the set of forests of height less or equal than n into subsets 
according to the number of vertices at each level, that is, 

H3Kx) := E \^\ xP = ^'W + E l-Fpl xP 
pev* P ev„ 

where we write x p as a shorthand for x^x^ 1 .... 

We let dp = d^^dxl ■ ■ ., and we consider the mapping B~ l from V into itself denned for 
any k > by 

B- 1 : pe V k ^B-\p) = (l,p) = (l,po,pi,...) G Vjfc+i 

Since forests f of height are characterized by the number of roots in f , we clearly have the 
formula 



(VpeVo |^p| = l) 



(-o) = r ^ = E^o 



p>0 



Proposition 3.11. 

Vn 



^ = n (iz^sO^W'M* II (tz^w) 



SpH^-^O) 



pev*_! 



pev n _i 



In addition, the generating function associated with the number of forests with multidegree 
p G V n; and 6mZi wi/t trees wit/j height n is given by 



Vn>i E K^))n^ P |xP= n (tt^t 

pev n pev n _i ^ x 



(p) 



l (0) 



Before getting into the proof of this result, we notice that the formula stated in the above 
proposition makes the Hilbert series computable at any finite vertex order, and at any height 
using any formal computation software. The first two orders can be handly computed. For 
n = 0, we have already seen that 

E£(x) = - J— = < => VP G V d p M°f(0) = \F P \ = 1 

X ° P>0 



By the recursion formula we find that 

hWx) = — *— x TT — - — 



n e ^ x ™) fcn 

n>0 \fc„>0 



(fen)n>0 

= E x o° ~*~ E i^pi ^o ^! 1 



po>o peVi 
This readily yields that for any p G Vi, we have 



I^pI = #< 



(k n )n>o G N N : E fc n = Po and E nk n =pi 
n>0 n>0 



31 



More generally, the (n + l)-th order Hilbert series are given by the formula 

= n (t^W" pI 

pev* v 

Recalling that for any m > 1, we have the following formal series expansion 

1 V m ( m -i + fc)! ^ 



1 — u ) ^-^ (m — 1)! fe! 



we find that 



P*+^M = TT V (l-^pl - l + fcp)! fcp s-i(p) 

11 ^ — i)i fc p ! 

pev*fe p >o Vl pl ; p 

l^p[-l + fc(p))! \ vEo6Vi fc(p) B -i(p) 

(fol-l) 

This implies that that for any p G V n+ i, we have 



^ 11 (I^pI -i)! fc(p)! ' 



\v\ V" / tt (|^ q |-l + fc(q))! \ rqfi , 

1 pl = ^ 11 (i^i - 1)! m ( } 

£ qe v**(q) s-Hq)= P 
Now, we come to the proof of proposition 13. 1 ll 

Proof of proposition [5TTT1 We consider the set S{n) of all trees with height less or equal 
to n 

S(n) = T*(n) := U^< k < n T{k) = {0} U [U peV * :po=1 T p ) 

= {0} U (Uo<fc<n U pG v fc :p =l ^p) 

The set of forests generated by S(n) coincide with the set .F*(n) of all the forests with height 
less or equal to n. That is, we have that 

(S(n)) = F*{n) : = U_i< fc < n ^"(A;) = U_i< fc < n U peVfc F v 
To take the final step, we consider the mapping 

Vp G V„, Vf G T p x(f) = * v(f) = x P o°x P i 
From previous decompositions, and using the Hilbert series expansions ()3.5j) . we find that 

E x * (f) = E E i^i xP = E i^pI xP 

fe(S(n)) -i<k<n P ev k P ev* 

Recalling that 

B:peV^B(p)£Vn and |T p | = |S(T p )| = \F B(p) \ 
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we also have that 



n = — n (— T° 



P) 



i -x n n ( i - x p 

i-x n n (i-x b(x) b (p)) 

with the shift operator B((xk)k>o) = (^fc+i)fe>o- This implies that 

1 -i— r -i—r ( 1 \ '" Pp ' -i—r / 1 \ l-^i 



n i_ X f(f) n n (i-x^-hp)) n (i-xb^&o 

fes(n) o<Knpev t _! v 7 pev;_, v 

The end of the proof follows then from the Hilbert series formula (|3.5|) . 



p 



AT 
n.q 



Recall now that the signed measures in the Laurent expansion of the distributions 
are indexed by the coalescence degrees of forests. To compute or estimate these Laurent 
expansions, we are therefore interested in a more precise Hilbert series, namely the one 
taking into account, besides the multidegrees of forests, their coalescence numbers. 

Definition 3.12. We denote by ^ p [q] := T v H c _1 (q) the set of forests in J- p with a 
prescribed coalescence sequence q (we also use the convention Fo[0] = {$})■ 

We let Cjr(x, y) be the (multidegree) Hilbert series of forests associated to the partition 
of the set of forests of height less or equal than n, and prescribed coalescence sequences into 
subsets parametrized by multidegrees, that is 

C£(x,y):= E E l^ P [q]|x p y q 
pev* qec(p) 

In the above display, C(p) stands for the set of coalescence multidegrees 

Vp G V C(p) := c(^p) = {q G N N : and q < (B(p) - 1) + } 
with (B(p) - 1)+ = {{p k - l) + ) fc >i, for any p = (p k )k>o G V. 

Proposition 3.13. The multidegree Hilbert series of coalescent forests C^(x, y) satisfies 
the recursive formulae 

C^(x,y) = C^,y) [I II ( l x b-i (p)v J 

Using the same lines of arguments as before, this proposition can used to derive a recur- 
sive formula for the explicit combinatorial calculation of the number of forests prescribed 
heights, and coalescence multi-indices. Notice that the first two orders are given by 



Ci(x,y) = and Ct(x,y) = TT -, -y- 

l-x Q ^ V ' 11 n _ rnT Po 7I (po-i)+n 

u p >o v 1 x o x \ Vo ) 



Now, we come to the proof of proposition 13. 13l 



33 



Proof. We use the same notation as the one used in the proof of proposition 13.111 Firstly, 
we notice that T p = U qG c( p ) ^"pfq], and ht(q) < (ht(p) — 1), for any q £ C(p), and where 
ht(q) := inf {n > : sup m>ra q m = 0}, for any q = (q m )m>o- Also observe that 

^B-i(p) = ^i, P ) C(B-\p)) = c(T (liP) ) = { Po - 1} x C(p) 

We next consider the mapping % 

Vp £ V n Vf £ ^ p [q] X (f) = x^ f ) y c ( f ) = (« ...<") (y 90 yf . . . y*- 1 ) 

Using the Hilbert series expansions ()3.5j) . we readily check that 



n ; J, lYCil) = n E x * ,(g) y** ] = E ** ,S V<* 

feS(n) y f6S(n) \ge(f> 

Arguing as in the proof of pr op osition 13 . 1 1 1 we find that 



p y q 



E x^v(s) = i+Ei^i xP + EE E i^pMi* 

ge(5(n)> pev fc=i pev fc qG C(p) 

= E E i^pMI xP y q 

PGV* qeC(p) 

In much the same way, we have that 

n i _ x «(f) yC (f) = n n n ( r_ x i> v> 

feS(n) J u fc=ipev fc :p () =iqec(p) x ' y - 

Therefore, 

n 



feS(n) l_x"( f ) y c < f ) 

/ \ |^(i, S(p))[(Pl-l.-B(q))]| 

rifc=o n PG v fc :po=i n q ec( P ) ^.^p] ^pi-d J 



This yields that 

|^B-icp)[q]| 

feS(n) 17 A;=0peV fc _ 1 qgC(B- 1 (p)) 



n 1 _ x v(t) yC (f) n n n ^ ( x _ x b-^ P ) yq ) 

L (p) y q 



/ 1 \ l^ P [B(q)]| 

n n 

p6V*_ iqG C(B-l( P )) V 



The proposition follows. 



Notice that the Hilbert series technique can be developed to any order of refinement. 
For example, it could be used to take into account, besides the number of vertices or of 
coalescences at each level, the cardinals #(f), so that the coefficients of the Laurent expan- 
sion of the measures Qn,q could be read, in the end, entirely on the corresponding Hilbert 
series. We leave the task of expressing the recursive formula to the interested reader, and 
simply point out that the technique allows an easy, systematic, recursive, computation of 
the coefficients of the expansion of the measures Q^q a ^ an Y order, both in n and q. The 
observation can be usefull, especially in view of the systematic development of numerical 
schemes and numerical approximations based on Feynman-Kac particle models. 
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3.8 Some forests expansions 

In the present section, and the forthcoming one, we take advantage of the langage of trees 
and of the results obtained on their statistics to compute the first orders of the Laurent 
functional representation of Qn,g an d, respectively, to derive a natural generalization to 
Feynman-Kac particle models of the classical Wick product formula. 

Let us fix n and q > 4, so that the notation q denotes, once again, the constant sequence 
of lenght n associated to q. As we have already pointed out, there is only one forest in 
.Fq without coalescence, which is the product of q trivial trees of height n. There is only 
one forest in with only one coalescence, at level i, that will be written fi j. Its only non 
trivial tree is the tree with one coalescence at level % and two leaves at level n. There are two 
forests with coalescence degree 2 and the two coalescences at levels i, written f| j and ff j. 
The notation fg i denotes the forest with only one non trivial tree with coalescence degree 2, 
a vertex with coalescence order 2 at level i and its three leaves at level n. The notation 
denotes the forest with two non trivial trees with coalescence degree 1 and the coalescence 
at level i. There are four forests with coalescence degree 2 and the two coalescences at levels 
i < j, written f^ij, k = 1...4. The forest • has one non trivial tree with coalescences at 
levels i and j and its tree leaves at level n. The forest ff jj has one non trivial tree with 
coalescences at levels i and j and its three leaves at the levels j,n,n. The forest f|j„- has 
two non trivial trees with one coalescence at level i, resp. j and their two leaves at level n. 
The forest t^ij ^ as ^ wo non trivial trees with one coalescence at level i, resp. j and their 
two respective leaves at level j, n, resp. n, n. 

Expanding the formulae for d l Q qtn and using the formulae obtained in Thm l3~%l for the 

cardinals #(f) and using that s(q,q — 2) = f ^ J 22pL ( se e for instance j^j, on page 63), 
we get the following result. 



1 N are given by the following formulae 



Corollary 3.14. The first three order terms in the polynomial functional representation of 
/(<?-!) \- / A fi,» 



d 2 



0<k<n 



Zn,q 



4^ E {& + ! (» - 3) *St -!<.-D & + ^ 

>' ' 0<k<n ^ ' 

E (lff«-(^+Aj#)) 

' 0<k<l<n 

(q-1) (AjL^t /A f w , A f ^ , (9-2)(9-3) .f|, fe , 



£ (a£< ; + (g - 2) (A^ + AS M ) + (g 2) 2 (g 3) AS fc '') 

0<k<l<n ^ ' 

3.9 A Wick product formula on forests 

Let now B S Q m (En) C Bb(En) be the set of symmetric functions F on E n such that 
(Di q -i®Jn){F)(xi,...,x q -i)= / F(a?!, . . . ,x q -i,x n ) J n (dx n ) = 
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Notice that B^ m (En) contains the set of functions F = (P) sym , with P G Poly(En), 
where Poly (En) C B b (En) stands for the subset of polynomial functions of the form 

p = ^ c(a) f a with f a = (f a ^ <g> . . . <g> f a ^) 

ael 

In the above display, / is a finite subset of N^, c G M 7 , and the elementary functions 
/ p G Bb(E n ) are choosen such that ~f n (f p ) = 0. For instance, we can take 

/ p = (S p " ?MG? P )) with g p eB b (E n ) 

Let us recall that, when explicited in the particle models setting, the traditional Wick 
product formula reads, for any integer q and any F G B S Q m (El) 

&Qo, q (F) = 0, i < | 

and, if is even, 

where f is the forest in Tq^ containing | copies of the tree of unit height, with two vertices at 
level 1 , and | copies of the tree with the root as unique vertex. For symmetric tensor product 
functions F = (f 1 <8> . . . <g) f q ) sym , associated with a collection of functions f % G B b (Eo), such 
that rjo(f l ) = 7o(/ 1 ) = 0, for any 1 < i < q, we readily find that 

^ \{i,i}63: 9 / 

as soon as q is even. In the above displayed formula, X q ranges over all partitions of [q] into 
pairs. In a more probabilistic language, the above formula can be interpreted as the q-ih 
order central moment 

d q i®o, q (F)=nwo(f 1 )...w (n) = Y J ( n nw (f)Wo(f))\ 

i q \{i,j}ei q ) 

of a Gaussian field Wq on the Banach space of functions B^Eq), such that for any pair of 
functions (<p,tp) G B b (E ), E(W (<p)) = 0, and E(W ((p)Wo(ip)) = r] ((pip). In terms of the 
Laplace moment generating function 

x = - exp Qx'C (/)x) = ((27r)9| J Q(/)|)1/2 I ^exp (-^(/rV + y'x) 

with the gaussian covariance matrix Co(f) = (Co(f l , P))i<i,j<q, we also have that 

aq / 1 \ F)1 B B 

^Q.„(F) = ^ e xp(- X ' Co (/)x)^ with —^^j...— 

In the present section, we will show that the Wick formula generalizes to forests in 
J~ n ,q °f arbitrary height. Let us start by listing various straightforward properties of trees 
and forests. A tree T with coalescence degree d has exactly (d + 1) leaves. A forest with 
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coalescence degree d has at most d non trivial trees, and the equality holds if and only 
if all its non trivial trees have coalescence degree 1. In particular, if a forest in T nq has 
coalescence degree d, it has at most (2d) leaves belonging to non trivial trees so that, if 
d < |, there is at least one vertex at level n + 1 belonging to a trivial tree (that is, a tree 
with coalescence number 0). 

The same reasoning shows that, when d = |, a forest in !F n ,q with coalescence degree d, 
and coalescence sequence r, does not contain a trivial tree of height (n + 1) if and only if 
it is the forest f r := Tq°Uq° ...T^ n U^ n , where we write T k for the unique tree of coalescence 
degree 1 with a coalescence at level k and its two leaves at level (n + 1), and where we write 
U k for the trivial tree of height k. 

We conclude this series of remarks by noting that, if f G J~nq can be written as the 
product (or disjoint union) of a forest g in T n . q -\ with f7 n +i, the trivial tree of height 
(n + 1), then, for any F £ B s b ym (En) we have, by definition of the measures 

Ai iq (F) = Al q (D lq _ 1 ® ln )(F) 

which is equal to if F £ B s Q ym {El). 

We are now in position to derive the forest Wick formula. 

Theorem 3.15. For any even integer q < N and any symmetric function F £ B$ y (En), 
we have 

Vk<q/2 d k Q n>q (F) = and d^ 2 Q n , q (F) = £ A% q F (3.7) 

r<q,|r[=f 

For odd integers q < N, the partial measure valued derivatives d k are the null measure on 
B S Q m (En), up to any order k < [q/2\. 

We close this section with a gaussian field interpretation of the Wick formula (|3.7|) . We 
further assume that q is an even integer. We consider a collection of independent gaussian 
fields (W k ) < k < n on the Banach spaces (B b (E k )) < k < n , with for any (ip k , ip k ) £ B b (E k ) 2 , and 

< k < n 

E{W k (^ k )) = and E(W k (<p k )W k (iP k )) = j k (<pM 

We also introduce the centered gaussian field V n on Bb(E n ) defined for any tp n £ Bb(E n ) by 
the following formula 

V n (<p n ) = ^ VTfcOO W k (Qk,n(<Pn)) 
0<k<n 

Let (ipi)l<i<q e Bb(E n ) q be a collection of functions such that 7n( ( / 3 «) = 0, for any 

1 < i < q. For the tensor product function F = \ J2aeg q (W(i) ® • • • ® ( Pa(q))i one can check 
that 

A ^ (i?) = e ^(!) rfc n ik(Qk,n(<Pi)QkAw)) 

ieio<k<nJ k £i k { (i,j)eJ k 

In the above displayed formula, the first sum is over the set X of all partitions I = (I k )o< k < n 
of [q] into (n + 1) blocks with cardinality \I k \ = (2r k ), and the second sum ranges over the 



37 



set Ik of all partitions of the sets into pairs, with < k < n. By definition of the 
gaussian fields (Wfe)o<fe<ni an d due to the classical Wick formula, we find that 



^A% q (F n ) = e n e ( n w k (QkA<Pi)) 



IelO<k<n \iei k 



= E E II IIv^^WMfe)) 

\Iel 0<k<ni&I k 

from which we arrive at the following formula 

E j^ri A nA F n) = E ( e e n n^j^M^)) 

r<q,|r|=| ' \r:|2r|=g / 0<k<niel k ) 

Recalling that all gaussian fields (Wfc)o<fc<n are independent and centered, we prove that 

d"/ 2 Q n , g (F n )=El n ( y, \/#)^(ft,n(<» I = e ( n k«) 

\l<i<<? \0<fc<n / / \1<«<9 

Written in a more synthetic way, we have proved the following formula 

d c >/ 2 q n , q (F n ) = E(v^(F n )) 



This result can alternatively be derived combining the L^-mean error estimates presented in 
[1J Thm.7.4.2], with the multidimensional central limit theorems presented in jl| Prop. 9. 4.1]. 

l<i<q 

converge in law to (V^((/4))i<i<<?- By the continuous mapping theorem, combined with sim- 
ple uniform integrability arguments, one checks that 



More precisely, the q-dimensional particle random fields (V^ v (^)) 1<{< := ( VNj^ {ip l ni 



iW 2 ®% q (F n ) = E ((V n N r q (F n )) ^ d^ 2 Q n , q (F n ) = E {V^(F n )) 



4 Extension to path-space models 

In the present section, we extend our previous analysis to the statistical study of path spaces. 
Due to the mean-field definition of the trajectories, it will appear soon that distributions 
on the path space are parametrized by forests with a slightly more complex structure than 
the classical ones. Namely, we have to introduce colored trees and forests (black and white, 
with a particular structure deviced to reflect the geometry of paths). These new objects are 
introduced and studied in the first two subsections; applications to the path space and to 
propagation of chaos properties are postponed to the last three sections. 



4.1 Colored trees, forests, and jungles 

Definition 4.1. A colored tree is a tree with colored vertices, with two distinguished colors, 
say black and white. Only black vertices may have an ingoing edge. That is, equivalently, all 
white vertices are leaves -notice that the converse is not true in general. A colored forest is a 
multiset of colored trees. The sets of colored trees and colored forests are denoted respectively 
by T and T . 
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Figure 3: A colored forest 

Notice that a colored forest can also be viewed as a commutative monomial over the set 
of colored trees. In this interpretation, the generating series techniques that we have devel- 
opped to deal with the enumeration of forests will apply to colored forests. The computation 
of the corresponding series is left to the interested reader. 

Most of the notions associated to trees and forests go over in a straightforward way 
to colored forests and colored trees. As a general rule, we will write a line over symbols 
associated to colored trees and colored forests. For instance, we write v k (T) = (v k (T), v' k (T)) 
for the number v k (T), resp. v' k (T), of white, resp. black vertices in the colored tree T, at 
each level k > 0. 

We also let V be the set of all sequences of pair of integers p = (p k ) k >o £ (N 2 ) N , with 
Pk = (Pk,p'k) f° r every k > 0, satisfying the following property 

3ht(p) G N s.t. inf p k + p' k > and sup pk + p' k = 

k<ht(p) k >ht(p) 

For any n G N, we denote by V„ C V the subset of sequences p such that the height of p, 
ht(p), is equal to n. Finally, we let Tp be the set of colored trees T with vu{T) = p k white 
vertices, and v' k (T) = p' k black vertices, at each level k > 0. Since T is a colored tree, this 
implies that po = and p' Q = 1, excepted if ht(p) = 0. In that case, po may also be equal 
to 1 (and then p' = 0) . Notice that the set of black and white vertices at each level k is in 
bijection with the disjoint union 

[Pk} = KPk,p' k )} :=bk]lM 

of the sets [p^] and \p' k }. 

For p G V n+ i, let us call, by analogy with the uncolored case, planar colored trees T of 
type p any sequence of maps (ao, . . . ,a n ), where 

a k = Hua'fc) G [p' k ] [pk+ Mp' k ] [P ' k+l] =■■ \p' k ]toWk + iK 

is the disjoint union of the pair of weakly increasing maps a k € \p' k ]^ Ph+1 \ and a k G bfc]^ fc+1 ^ 
In the above definition, the parameter level k runs from 0, to n, and we have used the 
conventions (0u0) = 0, and (au0) = a. Writing VT for the set of planar colored trees and 
VTp for the set of planar colored trees of type p, from the previous discussion, for any 
p G V ra+ i with po = 0, p'q = 1, we have that 

~ {(a , . . . ,a n ) G [lpi] x [p'jfe] x ... x K] [p " +l] : "fc = {a k ua' k ) with a k | and a' k |} 

A planar colored forest is, by definition, a sequence (or a non commutative monomial) 
of planar colored trees. 
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Since to any sequence {a k , a k ), a\ = (aj,[]a|,') G [p*,']K p M-i>Pfc+i)l i s associated a 

n n t 

unique element {a\ \J ■ ■■ \J a k > a k II ••• II a k') °f \P k ] LI ■■■ W[P k ] 1=1 ^ i=1 ^ ; planar 
colored forests are in bijection with sequences of disjoint union of weakly increasing maps, 
in the sense that (with self-explaining notations) 

VTp ~ lAp := {{a , ...,a n )eAp : a k = (a k ua k ) with a k | and a' k ]} 

for any p G V n+ i, with 

^p:=([Po] [Pl] xbi] fel x---x[K] [p " +l1 ) 

Definition 4.2. VFe associate with a sequence of mappings a = (ao, ...,a n ) G »4p, mi/t 
p G V n+ i, a colored forest F(a) of height (n + 1). T/ie vertex set of the colored forest 
F(a) is given by the disjoint union set (jJ^oiPkfj U (LRiotPfe]) ■ ^ e ec ^5 es correspond to 
the pairs (i,a k (i)), and (i',a' k (i')), with the pair of white and black colored vertices (i,i') G 
([pA;+i] U[p'a;_|_i]) • The sequence can be represented graphically by an entangled planar graph 
with colored vertices, or colored jungle , written J (a). Notice that edges between two levels 
are allowed to cross. At a given level k, all white vertices are drawn left to the black vertices. 
The set of colored jungles is written J . 

o • • 




o o • • • 



Figure 3: A colored jungle 

Once again, the notation we use on colored forests extends to colored jungles in a self- 
explanatory way. For instance, Jp, p G V Tt+ i, stands for the set of colored jungles j with 
v k(j) = Pk white vertices, and v' k (j) = p' k black vertices, at each level < k < (n + 1). For 
any p G V ra+ i, we also have that Jp ~ Ap. 

Our next objective is to make more precise the algebraic picture of these objects. For 
any integers p,p', q, q' , r, s, and any disjoint maps 

a = {am') G [r][( p ' p ')l and b = (b U b') G 

we let (a V b) be the disjoint union of maps 

(a u a r ) V {b U b') = {aVb) u {a V b') G [r + s ]^ +p '' q+q '^ 

Notice that the property of being weakly increasing is preserved by this construction. More 
generally, for any p G V, and q G V, and any sequences a = {a k ) k >Q G Ap, and b = 
(b~k)k>o G Aj, we set 

a V b = {a k V b k ) k > G -4p +q - (4.1) 

In the above displayed formula, we have used the conventions p + q = {p k + q k ,p' k + Q k ) k >o, 
a k = 0, and \ = 0, for any k > ht(p), and / > ht(q). Notice that, if a G TAp- and b G TA^, 
then a V b G TAp+q. 



40 



The main advantage of this algebraic structure lies in the following non commutative 
composition formula in the set of colored jungles 

7(5) 7(b) = 7(5 V b) G Jp+q 

Notice that any pair of permutations (s, s') G (<5 P x £/ p /) can be seen as a pair of mappings 
(s,s') G ^[p+p']^ x [p + p']^^- In this interpretation, the disjoint union map 

s = (s U s') G [p + P f p ' p 'V 

can be seen as a permutation map of the set [(p,p')] = ([p] Lib']) that preserves the sets [p], 
and [p']. The reverse assertion is also true. The set of these permutations is denoted by 

Q{p,p') '■= <3p\l<3p' 

We end this section with a technical lemma that permits to extend in a natural way the 
jungle automorphism group analysis developed in section 01 to their colored analogue. 

Lemma 4.3. For any collection of integers p,p' ,q,q' ,r,r' ,s, s' , for and sequence of permu- 
tation mappings ~p G G(p, p >), 9 G G(q, q '), ~<J G G( r ,r')> t £ G( s ,s')> an d f or an U V a ^ r °f disjoint 
mappings a G [r'][( p ' p '^, and b G [s'jtfe?)] we have the composition formula 

{a V T ) o (a V b) o (p V 9) = (a o a o p) V (r o b o 9) 

In addition, for any A G Q (p+q),( P ' + q ') we have that 

(a V b) o A G ([r'}^ V [ S '][^')l) =^ A G {G M V ^ } ) 

The two assertions are a direct consequence of the definition of the V product. In the 
first equation above, it should be understood that a G [r'][( p,p )l is viewed as a map in 
o ^ ( M U [ r '] ) " through the embedding of [r'\ into [r] ]J [r'] . The same observation will 
be used without further notice in the forthcoming paragraph. 

4.2 Automorphism groups on colored jungles 

As for the expansion of measures in section we will need a definition of equivalence classes 
of colored jungles, a parametrization of the classes and a computation of their cardinals. 
Thanks to the algebraic model of colored trees designed in section H~T1 most of the construc- 
tions in section 12 and section 01 can be extended without further work to colored jungles. 
Since all constructions can be mimicked, we only outline the main ideas of the generalization 
from jungles to colored jungles. 

We let p G V n+ \ be a sequence of (n+1) integer pairs p k = (pfc,p^,), with < k < (n+l). 
We associate with p the product permutation group 

Gp = Qp x g Vl x . . . x g Vn+i 

This group acts naturally on sequences of maps a = (ao,ai, a n ) G ^4p, and on jungles 
J(a) in by permutation of the colored vertices at each level. More formally, for any 
s = (so, G Gp this pair of actions is given by 

s(a) := (sooosf 1 , s^s"^) and sJ(a) := J(s(a)) 
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Two colored jungles in the same orbit under the action of Qp have the same underlying 
colored forest. Conversely, if two colored jungles in Jt> have the same underlying colored 
forests, they differ only by a permutation of the vertices of their colored graphs that preserves 
the colors of the vertices, and therefore are in the same orbit under the action of Qp. In 
other terms, we have the following lemma. 

Lemma 4.4. Equivalence classes of colored jungles in Jp under the action of the permu- 
tation groups Qp are in bijection with colored forests in Tp 

Notice the following corollary of the lemma: if f G Tp, with p G V n+ i, is a colored 
forest, then we can define unambiguously 

|f| = (|f|*)fc>o = (|a |,...,|a n |,(0,0),(0,0),...) 

for any choice (ao, ...,a n ) of a representative of f in J p. 

As for usual trees and forests, we write B for the map from colored trees to colored 
forests defined by removing the root of a colored tree, and we write also, as in section . 
B for the induced map from the set of colored forests into itself. 

Let us conclude this section by enumerating the number #(f) of colored jungles associ- 
ated to a given colored forest f G Tp, with p G V n+ \. The process is as in section and 
the result follows ultimately from the class formula 



#(?) 



p! 



\Stab{i)\ 

In the above displayed formula, p! stands for the multi-index factorial p! = (po\p' \...p n+ i\p' n+1 ty , 
and we have written abusively \Stab(f)\, for the cardinal of the stabilizer in Qp of any rep- 
resentative f' of f , where f is viewed as an equivalence class of colored jungles. 

Let us assume that f can written, as a monomial over the set of colored trees as 

i=l 

where the TjS are pairwise distincts and m\ > ... > m^. We write T\ for a set of represen- 
tatives of the Tj viewed as equivalence classes of colored jungles. 

As in section ETH we shall write S(B~ 1 (f)) for the unordered /c-uplet (mi, m^), called 
the symmetry multiset of the colored tree S _1 (f), and extend the notation to colored forests 
so that s(f) is the disjoint union of the S(Ti). The proof in section ETH goes then over without 
changes, excepted for the replacement of forests, planar forests and planar trees by their 
colored analogues. For instance, we have the recursion formula 

k 

\Stab(f)\ = Yl (mil \Stab B(Tl)\ m ^ 

8=1 

Expliciting the recursion in the Proposition, and following the wreath product analysis 
designed in section 131)1 we deduce the following theorem. 

Theorem 4.5. The number of colored jungles with underlying colored forest f G Tp, with 
p G V n+1 , is given by 

#(f) = - 



p! 



n s(B*(f))! 

;=-l 
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where we use the usual multi-index notation to define s(B l (f))\. In addition, we have the 
wreath product representation formula 

Stab{?) ~ g mi < Stab B^) > V . . . V Q mk < Stab B(T' k ) > 
4.3 Feynman-Kac semigroups 

In this section, the parameter n > represents a fixed time horizon. For any sequence of 
integers q = (qo, ... ,q n ) € N n+1 , and any — 1 < m < n, we set 

q'm : = Qm+l + q' m +l = Qk (4.2) 

m<k<n 

Notice that q' m = \[q m+ i}u[q' m+ i]\. We associate with q the unnormalized Feynman-Kac 
measures 

^ = 7^ qo ^--^^GM(B^ where E^^x-x^ 

is equipped with the tensor product sigma field. Points in Ein are indexed as follows 

y\d, , . . . , x ),..., {x n , . . . , x n jj 

Definition 4.6. We let (7p )o<p<n> and (Qp)o<p<n, be the collection of measures and integral 
operators defined by 



7 c, = u^.®...®^)®^*-^) (= 7 ^,., qp - M ' P r 



Q P q = {D lvo ® ■ ■ ■ ® D lqp _ 2 ) ® Q^- 1 '^ 

with the operators Q p 9p 1,<?p_1 - > from £^L 1 1+9p_1 into E q JS\ x E^' 1 defined by the tensor 
product formula 

Notice that, for any < p < n, jp is a positive measure on the product space 

E^^.E^x-.-xE^xE^ 
and Qp is a positive integral operator from E^_ x into Ep. 

Lemma 4.7. For any sequence of integers q 6 N n+1 , and any < p < n, we have 

In Ip ^*p,n 

with the semigroup (Qp 1 ,p 2 )i<p 1 <p 2 <n defined by 

Qpi,P2 = Qpi+i (Qpi+i^) ( = Qpi+i • • • QP2-1QP2) 
cmo? t/ie convention Qpi, Pl = Id, the identity operator, for p\ =P2- 
Proof. We start with observing that 

= 7o 090 ® 7? 91 ® ■ ■ ■ ® 7^2" 2 ® (7?V + *' ) ® Q®**)) 

Therefore, we find that 7n = 7^'_ 1 Qn- Using a simple induction, the lemma follows. ■ 
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4.4 Unnormalized particle measures 

In this section, we derive a functional representation and a Laurent expansion of particle 
tensor product measures on the path space similar to the ones obtained in Section 2.2. 

In the further development of this section, the time horizon n > is a fixed parameter, 
and we let q = (go, • • • , Qn) £ N ra+1 be a given sequence of (n + 1) integers. Before we turn 
to the study of the particle tensor product measures 

1% N = (7o) m ° ® (7i ) mi ® • • • ® (7^)® 9 " 
and the associated path space measures Q^ q G TW(E^) defined by 

< q : F G fi b (E2) —> Q^ q (F) = E {^ N (F)) G R 

let us introduce some useful definitions. We associate with q the pair of integer sequences 
q, and q' defined by 

q : = (lk,q'k)-i<k<n and q' := (q' k )-i<k<n 

with the convention g_i = 0, and the integer sequence q' introduced in (|4.2j) . Notice that 
Qn = °; 111 = 9o + 9o' ^(q) = ( n + and 

-l<fc<n 0<fc<n 0<k<n 

Definition 4.8. We let A nq 6e £/ie nonnegative measure valued functional on colored jungles 
J (a) G jTq, mi/i a = (ao,oi, ...,a„) G .4 q; and defined by 

A n , q : 7(a) — > A» q := i^D^CtfD^ . . . Q^D^ G A<(E«) 

uni/i i/ie Markov transitions Djjf- , resp. -Da fc , witha^ = {ak\ia' k ) G fefeluta^l^ from E^, 

resp. E 9 k k+Qk , into itself, and defined by 

D?, = Di <8...<8-Di <8Au. 
and /or any function F G B b (E q k k+qk ) by 

D- ak (F)((xl . . . ,xf ), Ofii, • • . ,yf)) = F((s?« . . . ( y < (1) , . . . 

Definition 4.9. We write B^ ym (En) for the elements F of Z3j,(E£J) i/iai are symmetric in 
the first qo variables, the next q\ variables,..., and the last q n variables. That is we have that 

V(<t , • • • , o-n) € [Q qo x ... x g qn ) F = (D ao <8 ... <8 D„ n ) F 

We consider the subset BQ Vm (En) C B^ ym (En) of all functions F such that 

(L> l!o ® . . . ® L> 1?p i (i^ 7p ) ® D lgp+i ® . . . (8 L>i 9n ) F = 

/or any < p < n. 
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Notice that £>g ym (E^) contains the tensor product 

(B s ym (E«°) ® . . . ® B s Q ym (E^)) = {F ® . . . ® F n : VO < p < n F p G B?^(^)} 

of the sets B l Q )m (Ep p ) introduced in the beginning of section with < p < n. For the 
same reasons as in equation 11.151 we also have the symmetry invariance property 

VFe£f m (E£) VsGSq VfGJq A* q (F) = A«g(F) 

By construction, we also have that 

VF G By m (E£) A^ q (F) = 

as soon as the colored forest f contains at least one trivial colored tree with a white leaf. 

Recall that a tree, or a colored tree is said to be trivial if its coalescence sequence is the 
null sequence of integers. Notice also that for any 1 < q', the set of colored forests with 
exactly lk(& [0><7fc_iD coalescent edges at level k, is given by 

Fq[l] := {fGfq : |f | = q' - 1} 

Thus, for any r < q', the set of colored forests .Fq-(r) with less than coalescent edges at 
level k is given by 

Fq(r) := Ui< r ^q[l] 

The coalescence degree of a colored forest f G Fq[l], is the sum |q' — 1| of the coalescence 
orders of its vertices. 

Notice that a colored forest with coalescence degree d has at most d non trivial colored 
trees. In addition, a colored forest with a coalescence degree d, has at most (2d) leaves 
belonging to non trivial trees. Since a colored forest in Fq has exactly |q| white leaves, if 
|q| > (2d), then it contains at least one trivial colored tree with a white leaf. 

Next, we discuss the situation where |q| is an even integer, and we characterize the 
subset of forests (in Fq), with a coalescence degree d = |q|/2, without any trivial colored 
tree with a white leaf. This charaterization follows the same lines of arguments as the ones 
presented on page E3 For any integers < k < I < m < n, we let Tk,i, m be the unique 
colored tree, with a single coalescence at level k, and white leaves at the levels (I + 1), and 
(m + 1). Notice that colored forests with (2rk) white leaves, no black leaves, and coalescent 
degree r^, with different pairs of coalescent edges at level k are necessarily of the form 

f (tfc) T~T 7p tk,l,m 

k ~ 11 k,l,m 
k<l<m<n 

for some families of integers t k = (tk,i, m )k<l<m<n such that |t fc | = Yl tk,l,m = r k- Therefore, 

k,l,m 

a colored forest in Tq, with coalescence degree |q|/2, a coalescence sequence r such that 
|r| = |q|/2, without a trivial tree with a white leaf, has necessarily the following form 

e (t) _ fCtoJj-rro f( tn )rrr n 
1 — T U Q ■ ■ ■ 1 n u n 

for some sequences of integers t = (tfc)o<fc<n! such that |tfc| = r^, for any < k < n. In the 
above displayed formula, Uk denotes the unique trivial tree with a single black leaf at level 
k, with < k < n. We write t! for J J ifc,«,m- 

k<l<m<n 

We are now ready to extend the Laurent expansions, and the Wick formula presented 
in theorem 12.71 and theorem 13.151 to particle models in path spaces. 
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Theorem 4.10. For any q = (qo, . . . , q n ), with |q| < N , and any F £ Bt v (E^) ; we have 
the Laurent expansion 



l<fc<|(q'-l)+| 



(4.3) 



ura£/i i/ie signed measures d Q n „ defined by 



a'Q„„(F)= E E '""'t.' #(T) aUf) 



r<q', |r|=fc fe^(r) 



(q'),f 



In addition, for any even integer |q| < N, and any symmetric function F £ £>Q ym (E£[) ; we 
have 



VA:<|q|/2 d k Q n ^(F) = and d^ 2 Q n ^(F) = £ £ 



q! ?(*) 

1 rA^ „F 



r<q',|rj = |q|/2 ||t||=r 



2<5(t) t ! n 'i 



with the integer sequence ||t|| = (|t^ | )o<fc<n? and the sum of the diagonal terms 5{t) = 
So<fc<Kn^fc,U- For °dd integers |q| < N, the partial measure valued derivatives d k are the 
null measure on BQ Vm (En), up to any order k < [|q|/2j . 

Before getting into the details of the proof, we mention that the Wick formula stated 
above has a natural interpretation in terms of the gaussian fields (Vfc)o<fc< n introduced in 
section EH Following the discussion given in that section, we get that 

£> |q|/2 Q n , q (F) = E ((V<f 90 ® . . . ® K® 9 ") (F)) 

for any tensor product function F of the following form F = Fq^ (g) ... (g) F n ^ n , with 

Fk, qh = E„ k eg qk (<Pl kil) ® • • -® ^ fefe) ), ¥>1 e B 6 (^ fc ) ; and 7& (<4) = 0, for any 1 < i < ?fc , 
and any < k < n. 

Proof. By definition of the particle model, and arguing as in the proof of (|2.1j) we find that 



E 



([(T^i)^ 1 ®^)®*] (?) 



AN) 
Sn-1 



JV \g>(g n -i+?i0 



= (7n-l 

for any p G B s b ym (E^ 1 x £$ n ). This yields that 

E ^'V) ) = 7*20*D^ (F) = 7 ^Q^ LN (F) 

from which we readily conclude that 



D 



dp) 



E {ln' N (F)) = 1(7^5 
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A simple induction yields that 



jVlq'l E ( q /)' j A «,q( F ) jylq'l E #( f ) A n,q( i? ) 



from which we find the following formula 

l<P<q' KH ' lPl feFq:|f|=p 

Using the Stirling formula (|2..3[) . we readily check that 

C = E E ^m^i^ E #(f) A L 

l<l<q'l<p<q' VH; lPl feFq:|f|=p 

From previous computations, we conclude that 

C=Ei E KlfU'-r)^ #(f)A^ 

Finally, we notice that .Fq-(O) reduces to the single class of all sequences of bijections in Aq. 
The end the proof of the first assertion is now clear. To end the proof of the theorem, we 
notice that (with t as in the expansion of (?' q '/ 2 Q n q (i ? )) 



tk,l,m ] -) ( E[fc<Kn 2 tfe,! ' i 



#(f 



Sincewehaves(|f (t) |,q'-r) = 1, and (q'Uw, = U.-i<k<n(l'k)q' k -r k +v for any r = (r fe ) < fc <„, 
we conclude that 

s(|f (t) |,q'-r) #(f(t ),_ q! 



with t! = Ilo<KKm<n tk,l,m," The end of the proof of the theorem is now straightforward.i 



4.5 Propagations of chaos type expansions 

This section is essentially concerned with applications of the differential forest expansion 
machinery developed earlier to propagations of chaos properties of interacting particle mod- 
els. In order to state, and prove the main results of this section, we need to introduce some 
notation. We shall work throughout with a fixed time horizon n > 0, and a constant particle 
block size q < N. We let be the set of integer sequences p = (pk)o<k<n £ N n+1 such 

that |p| = q. We associate to a given p = (pk)o<k<n G N n+1 the pair of integer sequences 
p' and p + q defined by 

V0 < k < n p' k := y~] pi and p + q : = (p , . . . ,p n -i,p n + q) 

k<l<n 
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We also denote by (G}-)o<k<n, an d G n p the collection of functions defined by 



Gk := -4tt x (vk(G k ) - G k ) and G® p := G o 0po 



lk(G k ) 

Finally, we consider the integral operators S'n P+ '^ from En +<? into E^ +i , defined for any 
function F n+ i £ Bj,(E^ +1 ) by the following formula 



S^(F n+1 ) := GT° 



_1 ® (^G„ ® 



sym 



with 



^ " 7n(l) 



- (F„ - ^(F n )) 



Proposition 4.11. For any q < N, and any n > 0, we have the polynomial decompositions 

.'-V ._ 17 / / i .A',/-, , ,. 



E» n := E ((1 - 7^(G„)/ 7 n(G n )) y j £ If d k E, 

q/2<k<(n+l)(q-l) 



q,n 



The derivatives of order q/2 < k < (n + l)(g — 1) are given by the following formula 

p GN £+lfc<l(p'-l) + l 

Proof. We first check the following decomposition 

l-7^(Gn)/7n(G„)= £ 7^) 

0<p<n 



(4.4) 



To prove this formula, we simple notice that 

7^(Gn)-7n(G n ) = 



0<p<n 0<p<n 

7^(G n - Vn (G n )) + [7^-i(G n -i) - 7«-i(G„_i)] x % (Gy 



0<p<n 



n %(ga 

p+l<fc<n 



with the convention JTj = 1. Finally, combining the multinomial decomposition 



(l-^{G n )h n {G n )) q = £ ^7^(GQ 



|P|=9 

with the Wick expansion stated in Thm l4.lT71 we conclude that 

E £C(^) = E E $ 

pGN™ +1 ^ peN^ +1 <?/2<fc<|(p'-l) + | 



(4.5) 



This completes the proof of the proposition. 
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Theorem 4.12. For any n > 0, the sequence of probability measures (JP^ +1 q )N>q i- s dif- 
ferentiable up to any order with d°¥ n+ i^ q = f?®+ 1; and the partial derivatives given by the 
following formula 

f) k w , - V ((g - 1) + |p|)! rsk^r „ p+q 
peu <« 2fc Nf+ 1 
At any order N > q, we have the exact formula 

l<k<[{N-q)/2\ 

Proof. We let F n+ i be a bounded measurable symmetric function on E^ +1 , such that 
^n+i (-fn+i ) = 0. By definition of the particle model, we have that 



E 



(tei) 09 ^!)!^) = ^(G„rx(^)^QSi(^i) 

= (l-O^X^OSl^n+l) 



with the sequence of random variables := (l — 7^ (G n )/ r y n (G n )^j . On the other hand, 
for any {!}, and m > 1, we have the more of less well known decomposition 



v ; 0<fc<m l<k<q+l 

with (g + = {q + k)\/q\. The proof of this formula is essentially based on the fact that 
for any n > 0, and u 7^ 1, we have 

/(„) = _L_ = y u k + and^ = n! 

JK 1 l-u ^ 1-u du n J 

0<k<n 

This implies that for any m > n, we have 
d n f n\ 



V- / , \ fc _ n d n (u m+1 \ 



<9u™ (1 - u) n+1 

n<k<m 

Recalling the Leibniz binomial derivation formula 

9™ , . , f n\ d n - k , „ N a fc 



k ) d^ {f) ft?**) 



dw 

0<k<n 



we find that 

1 d n fu m+l 



v r n ^ fm + 1), u( m+i )- fe (n ~ /c)! — 

n\du n \l-u) L> \k) {m+L)kU n! (1 -«)(»-*)+! 



0<fe<rt 



k J (l_ n )(n-fe)+l 
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and 



Thus, 



1 d n fu m+l 



n! du n V 1 — u 



u 



E 

0<fc<n 



m+1 \ u k+l 
n-k ) (1 - u) k+1 



1 _ {k + n)\ u k 

(1 - ~~ ^ n! fc! 



" +u m - n y 



0<k<n 



u 



k+1 



m + 1 



n-k J (l-u) k+1 



We prove the desired decomposition with replacing m by (m + q), and n by q. Using the 
above decomposition we find that 



0<k<m 

with the remainder term 

Rm N n(F n +l) 

We also have from proposition 14. Ill 



1 _ 7^(G„) 

7n(Gn) 



g + m \ / 7 iy( G ) 
<<? \ k + m j \ yn(G n ) 



A - 1 

7n(G„) 



sup V^E(| 7 4 v (G n )- 7n (G n ) 

JV>m+l 



m+l\ m+1 



< OO 



Finally, using the multinomial decomposition (|4.5|) . we conclude that 



0<fc<m V 1 K ' 7 

E ((?-!) + *)* E ^rt +q) ' N (s^(F n+1 ) 



P eN" +1 



0<fc<m 

This yields, for any m + q < N , the functional expansion 



V- ((g-l) + |p|)' piV ,o' 

^ (o-l)! p! ^™>P+<r" + • 



j>7V 

n+1,9 

|p|<m 



< oo. The last asser- 



tv 



with a remainder measure n , such that supjv>i 

iV(m+l)/2 

tion follows from the regularity hypothesis II . II on the potential functions. This implies that 
for any k < (m + l)/2 we have 

|p|<m 
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Notice that fc-th order derivative measure d k Q n ^ p+q only involves colored forests with less 
than k coalescent branches, from the original root, up to the final level. If |p| > (2k), then 
these colored forests contain at least one elementary tree with a white leaf. By definition of 
the operator Sn +q we find that 

d k Q n>P+g ^ +q (F n+1 ) = 

This yields that 

P n+l,q = Vn+1 + X] Jfk ^n+l.g + ^n+l,q 

l<fc<((7V-g)+l)/2 

with a remainder measure q such that supjy >1 

of the proof of the second assertion of the theorem using the fact that 

[(N - q)/2\ <(N- q)/2 < ({N - q) + l)/2 

The Theorem follows. ■ 



We end this section, with a series of some direct consequences of the above theorem: 

• The differential forest expansions presented in theorem 14. 1 21 allow us to deduce precise 
strong propagation of chaos estimates. For instance, for any N > (q + 7) (so that 
(l(N - q)/2\ - 1) > 2) we have 



sup N 

N>q+7 



n+l,q Vn+l jy ¥n + 1 ' 



< OO 



with a first order partial derivative given by the formula 

d 1 F n+ltq (F n+1 ) = d 1 q n , q Q%(F n+1 ) + q £ d l Q nMm S^(F n+1 ) 



0<m<n 



with the sequence of integers q m = (l m (k) + ql n (k))o<k<n- The first term in the 
right hand side in the above displayed formula has been treated in corollary 13.141 on 
page 1351 and we have that 



q(q - 1) 



n+1, 

0<k<n 



jt 1] (d(x 2 , . . . , x*)) Qt q n+l {F n+1 ){x\ x 2 , x 3 , . . . , x«) 



0<fc<n 



Each of the terms 9 1 Q n qm 5'n m (i ? n+i) only involves the colored forests 



VO < k < m f fej7n := U k T k>m>n U^ +1 
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associated with the trees T^ mn , and XJ\, introduced on page 1451 After some elementary 
manipulations, we find that 

^ 1 Qn,q m Sn m (Fn+l) 

= q Y, &n,CS% m (F n+1 ) 
0<k<m 

= q E ^(1) [ ^\d(x\...,x^))Q k , m (G m )(x 1 )Qtl +1 (F n+1 )(x 1 



1 ■ ■ ■ l M 



A Wick formula derives from theorem I4.1UI More precisely, for any F £ B S Q Ui {E'/ l+1 ), 
the partial derivatives d k W n+ i A (F) are null up to any order k < [q/2\, and for any 
even integer q, we have 

d«/ 2 F n+1 , q (F) =7n+i(l)- g d^ 2 q n , q QT +1 {F) 

as soon as iV > 2{q + 2) (so that ([(N — q)/2\ — 1) > q/2). To prove this claim, we 
notice that for any k < q/2, and any F £ B^ /m (E^ l+1 ), we have 

Vp 6 U </<2feNf +1 k < q/2 < (|p| + q)/2 and therefore d k Q njP+q S% +q {F) = 

This yields that d k F n+ \ jq (F) = 0, for any even integer q, and any k < q/2. In the 
case k = q/2, we have k = q/2 = (|p| + q)/2 if, and only if, p coincide with the null 
sequence 0. 

We let ||//||,g=y» = sup FeB sym be the Zolotarev seminorm on M.{E^ l+l ) associated 

with the collection of functions 

^ : ={ feB nO ■ ll^ll < 1} 
For any even integer q such that (q + 2) < N/2, we have 

1 1 3i/ 2 n o® q 



sup N l+q l 2 

N>q+7 



7,N 
n+1, q 



< OO 

°0,1 



Combining the Wick formula stated above with the Borel-Cantelli lemma, we obtain 
for all q > 4 the almost sure convergence result 

lim ( V %) Gq (F) = p.s. 

for any bounded symmetric function F £ BQ fra (E q ). This result is an extension of the 
law of large numbers for [/-statistics obtained by W. Hoeffding [7] for independent 
and identically distributed random variables to interacting particle models. 

We mention that the same lines of arguments used in the proof of theorem 14, 121 show 
that the sequence of probability measures 



■ F£W + i)^K +1 , q (F):=EU V » +1 )^(F) 
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is differentiable up to order [(N — q)/2\, with <9°P n+li(J = 7/®?^ and the partial deriva- 
tives given for any 1 < k < \_(N — q)/2\ by the following formula 

8 k F n+1>q (F) = £ k ® n+Up>q) (Gl s F) 

peu < !<2fc N; i+1 

In the same way, for any F £ B'q' 111 (E 9 n+1 ) , and N > 2(q + 2), the partial derivatives 
d k Q n +i,q(F) are null up to any order k < q/2, and for any even integer q, we have 
that 

d^ 2 ¥ n+1>q (F) = 7 n+i(l)- 9 d^ 2 Q n+1 , q (F) 

• Finally, the Wick formula stated above allows to deduce sharp L^-mean error bound. 
To see this claim, we simply observe that 

F = (f-Vn+l(f)r q with / E B b (E n+1 ) => F% +hq (F) = E((^ +1 )^(f- Vn+1 (f))^) 
and for any even integer q such that (q + 2) < iV/2, we have the inequality 

1 1 



sup N 1+q ' 2 

N>q+7 



iN 

n+l « 7 n+l(l) 9 



< oo 

oaym 
°0,1 
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